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ABSTRACT 

We present 107 maps of continuum emission at 350 fi m from Galactic molecular clumps. Observed 
sources were mainly selected from the Bolocam Galactic Plane Survey (BGPS) catalog, with 3 addi¬ 
tional maps covering star forming regions in the outer Galaxy. The higher resolution of the SHARC-II 
images ( 8 .5" beam) compared with the 1.1 mm images from BGPS (33" beam) allowed us to identify 
a large population of smaller substructures within the clumps. A catalog is presented for the 1386 
sources extracted from the 350 ftm maps. The color temperature distribution of clumps based on 
the two wavelengths has a median of 13.3 K and mean of 16.3 d= 0.4 K, assuming an opacity law 
index of 1.7. For the structures with good determination of color temperatures, the mean ratio of gas 
temperature, determined from NH 3 observations, to dust color temperature is 0.88 and the median 
ratio is 0.76. About half the clumps have more than two substructures and 22 clumps have more than 
10. The fraction of the mass in dense substructures seen at 350 ftm compared to the mass of their 
parental clump is ~0.19, and the surface densities of these substructures are, on average, 2.2 times 
those seen in the clumps identified at 1.1 mm. For a well-characterized sample, 88 structures (31%) 
exceed a surface density of 0.2 g cm -2 , and 18 ( 6 %) exceed 1.0 g cm -2 , thresholds for massive star 
formation suggested by theorists. 

Keywords: star: formation — submillimeter: ISM — Galaxy: structure — surveys — catalogs — dust, 
extinction 


1. INTRODUCTION 

Although high-mass stars play a key role in the struc¬ 
ture and evolution of our Galaxy, there are still many 
questions that need to be answered before we have a 
paradigm for their formation process (Stahler et ah, 2000; 
Bally & Zinnecker, 2005; Zinnecker & Yorke, 2007; Tan et 
ah, 2014). Adopting the nomenclature of Williams et al. 
(2000), star clusters form in dense clumps (M ~ 50 — 500 
M©, sizes ~ 0.3 — 3 pc) of giant molecular clouds while 
individual stars form in cores (M ~ 0.5 — 5 M©, sizes 
~ 0.03 — 0.2 pc). An important step in the observational 
study of massive star formation is the identification and 
characterization of clumps and their resolution into cores. 
Studies of compact structures in our Galaxy will give 
fundamental “ground truth” for the analysis of nearby 
galaxies where individual clumps and cores cannot be re¬ 
solved, and provide source information for observations 
of the Milky Way with a new generation of instruments, 
such as the Atacama Large Millimeter/Submillimeter Ar¬ 
ray. 

There are complementary approaches to the study of 
massive star forming regions. We may consider obtaining 
molecular line observations, which are a good diagnostic 
of line-of-sight motions in a cloud, but variations in tracer 
abundances due to depletion and complex chemical pro¬ 
cessing, excitation conditions and the impacts of radia- 
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tion fields and shocks add uncertainties (Evans, 1999). 
We may also proceed by mapping the continuum emis¬ 
sion toward these sources. Newly formed stars heat the 
dust of their parental molecular cores and this dust emits 
in the far infrared and submillimeter wavelengths (Garay 
& Lizano, 1999). Continuum emission observations at 
(sub) millimeter wavelengths offer a reliable tracer of the 
column density and clump masses due to the low optical 
depth of the dust (e.g. Johnstone & Bally, 2006), and 
they have made an increasingly significant contribution 
to the understanding of star formation. Reddening and 
extinction of starlight provide also a measure of column 
densities (Bohlin et ah, 1978), and mid-IR large scale sur¬ 
veys have identified dark patches obscuring the bright- 
background, which are associated with dense molecular 
gas (e.g., Egan et al., 1998; Simon et ah, 2006). Stud¬ 
ies of the physical properties of these regions, called In¬ 
frared Dark Clouds, suggest they are the birthplaces of 
high-mass stars and stellar clusters (e.g., Rathborne et 
al., 2006, 2010; Battersby et al., 2010). In addition, near 
infrared extended emission has also been considered a 
tracer of the column density of dense dark clouds, inter¬ 
preted as scattered ambient starlight observed from the 
outer parts of dense cores (“cloudshine”, e.g. Foster & 
Goodman, 2006; Padoan et al., 2006). 

Several studies have made progress toward the char¬ 
acterization of massive star forming regions, but they 
are usually biased because target selection required ex- 
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isting signposts of newly form high-mass stars, including 
maser emission, far infrared emission, or the presence 
of H II regions (Churchwell et ah, 1990; Cesaroni et ah, 
1991; Plume et ah, 1992). Therefore, the identification of 
clumps with a less biased survey will improve our under¬ 
standing of the early stages of massive stars. Examples 
include large scale surveys of particular regions, such as 
the study of the W3 giant molecular cloud with Sub- 
millimetre Common-User Bolometer Array (SCUBA) 
by Moore et al. (2007) or the Cygnus-X molecular cloud 
complex region with MAMBO receiver by Motte et al. 
(2007). Large-scale blind surveys of the Galactic Plane 
at far-IR, submm and mm wavelengths have recently 
been completed. The ATLAS GAL survey mapped the 
Southern part of the plane at 870 fi m using the 12 m 
APEX telescope (Schuller et al., 2009; Contreras et al., 
2013). Space borne missions, such as WISE (Wright et 
al., 2010) and Spitzer j GLIMPSE (Benjamin et al., 2003) 
and MIPSGAL (Carey et al., 2009), have mapped most of 
the Galactic plane at a resolution of a few to 10", at wave¬ 
lengths between 3.6 and 24 /im. The Herschel/ Hi-GAL 
survey (Molinari et al., 2010) covered the entire Galactic 
plane with resolution between 5" and 36", in five bands 
at 70, 160, 250, 350 and 500 /im. In addition, the 15-m 
James Clerk Maxwell Telescope is using SCUBA-2 850 
/im and 450 /im emission to map the Galactic plane and 
several nearby Gould Belt clouds. These surveys have 
provided a robust tracer of dust properties (column den¬ 
sity and mean grain temperature), for different phases of 
the interstellar medium. 

The Bolocam Galactic Plane Survey (BGPS; Aguirre 
et al., 2011) is one of the first ground-based surveys to 
map the continuum emission at 1.1 mm over a large re¬ 
gion of the northern Galactic plane. The BGPS coverage 
region extended from l = —10° to +86.5° for b = ±0.5°, 
at a resolution of ^30", and it identified 8454 compact 
sources throughout the Galaxy (Rosolowsky et al., 2010). 
Depending on their distance, millimeter features may be 
cores, clumps, or clouds (McKee & Ostriker, 2007), al¬ 
though BGPS sources are mostly clumps (Dunham et al., 
2011). An improved reduction and a new version (V2.0) 
of the BGPS maps and catalog was released (Ginsburg et 
al., 2013), including additional ~20 square degree cover¬ 
age area in the third and fourth quadrant, and ^2 square 
degree in the 1st quadrant. This new reduction shows 
an improved flux calibration with respect to the former 
version (VI.0), in agreement with other data surveys on 
the Galaxy. The characterization of the angular trans¬ 
fer function of the Bolocam pipeline shows that the flux 
recovery is above 90% for scales between 33" and 80", 
with a steep drop for scales above ^100". While ver¬ 
sion 2.0 of the catalog contains 8559 sources, 35 sources 
were unintentionally excluded, and they were added in 
a new version (V2.1), for a total of 8594 sources. Fig¬ 
ure 1 shows the distribution of 40" aperture flux density 
of V2.1 sources. Although the basic statistical properties 
of the catalogs remains unaltered, there are some effects 
on the source extraction when considering the new map 
versions with improved quality. Comparing the source 
extraction in the V1.0 and V2.1 overlap region, V2.1 re¬ 
covers fewer compact sources (8004), with only ~ 70% of 
those having a clear VI.0 match, and these are consid¬ 
ered more reliable and stable identified sources (see Fig. 
10 in Ginsburg et al., 2013). Fig. 11 and 12 in Ginsburg 


et al. (2013) show visual examples of the lost sources, 
which are generally small sources on the shoulders of 
bright ones that have now been included in the bright 
source in the catalog. 

We selected a sample of sources from the BGPS cat¬ 
alog, initially from VI.0 and later from V2.1, towards 
the most crowded and brightest portions of the Galactic 
Plane including the Galactic Center, the Molecular Ring 
near l = 30° and Cygnus-X, and used the Submillimeter 
High Angular Resolution Camera II (SHARC-II; Dow¬ 
ell et al., 2003) to obtain continuum maps at 350 /im 
of them. The high resolution of the SHARC-II observa¬ 
tions (FWHM beam size of 8.5") is in particular very use¬ 
ful when we observe small structures at great distances 
in the Galaxy, and the comparison of flux densities of 
sources at 350 fi m and 1.1 mm will give us constraints 
on their dust temperatures, hence improving the deter¬ 
mination of masses and column densities of the molecular 
clumps. 

SHARC-II has been used before for studying low-mass 
dense cores from the Spitzer Legacy Program “From 
Molecular Cores to Planet-Formation Disks” (Wu et al., 
2007, A. Suresh et al. in preparation), concluding that 
observations with this instrument are better to distin¬ 
guish between starless and protostellar cores than ob¬ 
servations at longer wavelengths. While Mueller et al. 
(2002) studied dense clumps selected for having massive 
star formation using SHARC at 350 /im, there has not 
been a comprehensive study of BGPS-selected sources 
with SHARC-II. A recent study by Ragan et al. (2013) 
shows how ground-based observations at 350 /im contin¬ 
uum emission of regions of high-mass and cluster forma¬ 
tion can give information on substructures within them. 
For a set of 11 nearby Infrared-dark clouds (IRDCs), 
they obtained data on 350 /im continuum emission with 
the SABOCA instrument at the APEX telescope, finding 
a large population of small core-like structures. Higher 
resolution observations at submm wavelengths allow the 
study of the complex process of massive star formation 
from fragmentation of molecular clumps to individual 
substructures, such as cold and hot cores. 

We present a set of 107 continuum emission maps at 
350 /im obtained with SHARC-II. These include 104 
maps of clump-like BGPS V2.1 sources, and 3 maps 
of outer Galaxy sources not included in the catalog. 
This paper is organized as follows. Section 2 describes 
the selection of targets from the BGPS catalog and the 
SHARC-II instrument. Section 3 describes the data re¬ 
duction and calibration methods for our set of submm 
continuum maps. Section 4 presents the results of the 
extraction of sources in the images, including source re¬ 
covery tests, a catalog of these sources, the procedure for 
association between sources in the 350 /im maps and the 
1.1 mm sources, and the comparison between SHARC-II 
and Herschel at 350 /im. Section 5 describes the deter¬ 
mination of color temperatures for BGPS sources found 
in our set of maps,comparing them with temperature de¬ 
termination from spectral energy distribution (SED) fit¬ 
ting from Hi-GAL and other continuum surveys. Cor¬ 
relations between our dust temperature estimations and 
high-density molecular tracers are also presented. In Sec¬ 
tion 6, we explore consequences from our results for mass 
and surface density estimations of millimeter sources. In 
Section 7, we present conclusions from our analysis.. 
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2. OBSERVATIONS 

Observations were taken in several runs between June 
2006 and September 2012 at the Caltech Submillimeter 
Observatory (CSO 1 ). The observations required excel¬ 
lent weather conditions, with a measured optical depth 
at 225 GHz ranging between 0.030 and 0.094 (measured 
at zenith position), with an average value of T225 GHz — 
0.057, corresponding to an optical depth at 350 /im (857 
GHz) of T 350 /xm — 1-41. In general, we observed our 
targets when they were culminating, and we avoided ob¬ 
serving sources below an elevation of 30° or above 80°. 

2 . 1 . Target selection 

The targets were initially selected from sources in the 
BGPS VI .0 catalog. We chose sources with strong emis¬ 
sion at 1.1 mm (typically with peak emission above 1 
Jy beam -1 ) toward densely populated regions, many of 
them known as active star forming regions, such as the 
Galactic Center, l = 30°, Cygnus-X region, W3 and W5 
molecular clouds, GemOBl region, and others. Several of 
those targets were extended and they were likely to have 
unresolved substructure at a resolution of 30". Most of 
the time the observation fields were centered directly on 
a VI.0 source, but in some cases we tried to cover several 
sources in the same field. Starting in December 2009, we 
added to our list of targets bright sources selected from 
new Bolocam maps that were later incorporated in BGPS 
V2.0. Considering that almost all targets from VI.0 cat¬ 
alog are incorporated in V 2.1 (Ginsburg et ah, 2013), we 
consider the latest version of the catalog for our analysis. 
In just one of our maps, L136.52+1.24, there are three 
sources marginally detected in V 1.0 but not included in 
subsequent V2.0 or V 2.1 catalog versions. We also ob¬ 
tained 350 /im maps of the outer Galaxy star forming 
regions SH2-209, SH2-307 and IRAS 03233+5833 (e.g., 
Klein et ah, 2005; Richards et ah, 2012 ), which were not 
covered by the BGPS survey. The maps of these regions 
are L151.61-0.24, L234.57+0.82, and L142.01+1.77, re¬ 
spectively. These observations are included for complete¬ 
ness, but they are not included in analysis that requires 
BGPS data. 

For a better characterization of our targets with re¬ 
spect to the complete catalog, we consider the brightest 
V2.1 source in each map as the representative target for 
that map. In some cases, where adjacent maps share the 
brightest source, we consider the second brightest source 
as a representative target. From the 104 maps at 350 
/im with BGPS sources, we have then a sub-sample of 
V2.1 sources considered as representative targets. Fig¬ 
ure 1 shows in red the distribution of flux density in 40" 
apertures for our sample of 104 representative targets. 
The distribution has an average of 2.63 Jy and a median 
of 1.37 Jy, much higher than the values for the overall 
BGPS catalog. However, weaker sources were covered in 
the same fields, so the final flux density distribution of 
BGPS sources with 350 /im data is closer to that of the 
full catalog (see black line in Figure 1). We found a total 
of 619 V 2.1 sources in our sample of SHARC-II maps, 
with an average and median flux density of 0.79 Jy and 

1 The Caltech Submillimeter Observatory was operated by the 
California Institute of Technology, until 20l3 March 31 under co¬ 
operative agreement with the National Science Foundation (AST- 
0838261) 


0.29 Jy, respectively, about twice the values for the com¬ 
plete V 2.1 catalog (average of 0.35 Jy and median of 0.16 
Jy for 40" apertures). 

2 . 2 . Instrument and Observational Method 

SHARC-II is a background-limited facility camera op¬ 
erating at 350 /im or 450 /im, located at the CSO on 
Mauna Kea, Hawaii. It consists of a 32x12 array of 
doped silicon ‘pop-up’ bolometers and it has a mean 
beam size of 8.5" at 350 /im. The array has a full size 
of 2.59' x 0.97', which makes the instrument start to lose 
sensitivity to emission on scales larger than ~ 1 '. We used 
the box scan mode of SHARC-II without chopping 2 for 
mapping regions with sizes significantly larger than the 
size of the array. In the box scan mode, the scan pattern 
moves across the rectangular field of view, starting in the 
center of the array and going straight until the edge of the 
bounding box, where it changes trajectory emulating the 
bounce of a ball on a billiard table. The advantage of this 
scan mode is that it provides better isolation against 1/f 
noise, with more uniformly distributed noise across the 
field of view in comparison with the alternative sweep 
mode with Lissajous scan pattern. Different box sizes 
were used in our data acquisition, between 2.5'x 2.5' and 
lS'xlS', although the bulk of our observations were ob¬ 
tained with 10'xlO' boxes. Column 3 in Table 2 shows 
the size for each map. 

For all our observations we used the Dish Surface Op¬ 
timization System (DSOS ) 3 to correct the dish-surface 
figure for imperfections and gravitational deformations 
as the dish moves in elevation during observations. 

3. DATA REDUCTION AND CALIBRATION 
3.1. Data reduction 

The 350 fi m raw data were reduced using the Com¬ 
prehensive Reduction Utility for SHARC-II (CRUSH; 
Kovacs, 2006). CRUSH is a Java-based software which 
solves a series of iterative models, attempting to repro¬ 
duce atmospheric and instrumental effects on the ob¬ 
servations, removing the sky emission common to all 
bolometers at first order (e.g., Beelen et ah, 2006). We 
used version 2.10 of CRUSH, publicly accessible from the 
CSO webpage. 

Figure 2 shows an example of the error or noise maps 
that reduction with CRUSH produces for each observed 
region. In general, all SHARC-II maps are much noisier 
toward the edges by a factor of ^4 with respect to the 
measured average noise of the image. At a distance of 
^ 20 " (twice the beam size) from the edge, the noise is 
only increased a factor of ~2. Also, when taking multiple 
scans for a single map at different runs, the orientation 
of the field of view rotates, and the regions that do not 
overlap in the integration of the source show a significant 
increase in the noise. The rest of the map shows some 
residuals of the pattern used in the box scan mode, but in 
general the inner parts of these maps are well represented 
by uniform noise. 

3.2. Calibration 

2 See http://www.submm.caltech.edu/~sharc/ 

3 See http://www.cso.caltech.edu/dsos/DSOS_MLeong.html 
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The maps were reduced in instrumental units of /iV, 
and we estimated the flux conversion factors Cq in or¬ 
der to measure the flux densities of the sources in units 
of Jy. On our analysis, we followed the same method 
used on SHARC-II maps by Wu et al. (2007), and on 
SCUBA maps by Shirley et al. (2000), based on aperture 
calibrations over a given angular size 0. Flux estima¬ 
tion obtained from aperture corrected photometry are 
less sensitive to uncertainties due to effects of sidelobes. 
The observed flux density for an aperture of diameter 0 
is S u (A, 6) = Cq U(A, 0), with U(A, 6) the voltage 
measured at wavelength A in an aperture of diameter 0. 
We used the fluxes from Mars, Uranus and Neptune to 
obtain the flux conversion factors for two different aper¬ 
tures, 20" and 40". We restricted our flux estimation 
to 40" considering the limited sensitivity of the instru¬ 
ment to large scale emission. The bulk of the compact 
structures identified on the maps have sizes below this 
limit. 

To compare the 350 pm data with the 1.1 mm im¬ 
ages from BGPS survey, we also convolved the 350 p m 

to match the 33" FWHM beam from Bolocam. There¬ 
to// 

fore, we also estimated the flux conversion factor C 4 q 
for an aperture of 40" using the images of the calibration 
sources convolved to a 33" beam. The flux densities of 
sources and calibrators in instrumental units were mea¬ 
sured using Starlink’s aperadd package. 

In addition, we estimated the flux conversion factor 
for one beam, Cbeam, which we used to express the peak 
flux of each source and the 1 a rms of the maps in units 
Jy beam -1 . Cbeam is obtained from the integrated flux 
density of the calibrators in one beam and the peak pixel 
of the map in instrument units. We have then an aver¬ 
age value for the conversion pV to Jy pixel -1 = 32.23 
Jy beam -1 , for the beam size and pixel scale of the in¬ 
strument. Therefore, the value of Cq will depend of the 
number of pixels considered. From the maps, the average 
value of Cbeam = 8.27 Jy beam -1 /xV -1 , and then 

1 Jy beam -1 =503.8 MJy sr -1 
= 0.031 Jy pixel -1 
= 0.121 /iV -1 . 

For the calibration of the images, only the observa¬ 
tions of planets were considered. While the secondary 
calibrators mentioned in section 4.1 are useful to check 
the pointing when planets are not available, they are not 
as bright and their flux densities are not as well known as 
planets. Therefore, we avoided using them to minimize 
uncertainties in the estimation of calibration factors. 

Table 6 in the Appendix section gives the observation 
dates (column 1), the planet observed (column 2), the ob¬ 
served opacity at 225 GHz at that date (column 3), the 
computed values of Cbeam (column 4) and Cq at different 
apertures (columns 5 and 6). Column 7 gives the calibra¬ 
tion factor for an aperture of 40" obtained for the planet 
images convolved to a beam of 33". Table 1 shows the 
average values of Cq for the different observation epochs. 
According to these results, the uncertainties in the cali¬ 
bration of the maps are between 15-20%, consistent with 
previous observations and calibration estimations for the 
SHARC-II instrument (Wu et al., 2007). 

While some sources may have significant contributions 
to the broad-band flux density from line emission, it is 


usually a small fraction of the dust continuum emission 
(e.g., Groesbeck 1995). The continuum flux measured 
with SHARC-II (780 - 910 GHz for the 350 pm filter) 
could in principle be contaminated by the CO(7 —6) line. 
Hatchell & Dunham (2009) estimated a contribution less 
than 100 mJy beam -1 in the SHARC-II 350 pm band for 
a CO(7 — 6) integrated intensity of 100 K km s -1 , which 
is smaller that the typical sensitivity found in our maps, 
and they considered that the CO line may contribute 
^20% of the continuum emission in the 350 pm band. 
We adopted that limit in this work as a conservative 
value for possible line contamination. 

4. RESULTS 
4.1. Map description 

Figure 3 shows as an example the 350 pm contin¬ 
uum map of L133.71+1.21, corresponding to the W3 
Main region (see Rivera-Ingraham et al., 2013; Megeath 
et al., 2008, and references therein). Besides the two 
bright infrared sources IRS5 and IRS4 (Wynn-Williams 
et al., 1972), the small beam size of the 350 pm maps al¬ 
lows us to identify several sources near them, and some 
elongated, filament ary-like structures and scattered faint 
sources. These types of features are found recurrently in 
the 350 pm maps and we characterize their correspon¬ 
dence with 1.1 mm sources later in section 4.5. 

We obtained a total sample of 107 maps, with an aver¬ 
age noise of 744 mJy beam -1 , and a standard deviation 
of 443 mJy beam -1 . These maps are presented in Fig¬ 
ure 25 in the Appendix. Because the noise increases at 
the edges of the images (see Section 3.1), we estimated a 
more representative value of the mean rms of each map 
on a region toward the central area. In the 10' x 10' maps, 
the considered area has a diameter of 7'. Column 7 in 
Table 2 shows the representative rms noise of each map. 
The average la noise value estimated this way is 478 
mJy beam -1 , with a median of 375 mJy beam -1 and a 
standard deviation of 298 mJy beam -1 . Each map has 
an angular scale of 1.618 arc-seconds per pixel. Table 2 
lists, for each map observed, the configuration used, the 
size and center of each map, the date on which the map 
was taken, and the 1 a noise of the map in units of mJy 
beam -1 . Integration time was obtained in blocks of ~ 
14 minutes. The pointing was checked on planets such 
as Mars, Uranus and Neptune, and with secondary ob¬ 
jects when the planets were not available, such as IRAS 
16293-2422, K350, G34.3, W75N and CRL 618. The 
blind pointing uncertainty varied between 1.8" and 2.1" 
for azimuth and between 0.8" and 1.6" for zenith angle. 

4.2. Source Extraction 

To facilitate comparison with the BGPS sources, we 
used the same algorithm for source extraction as was 
used for the BGPS catalogs, Bolocat (Rosolowsky et al., 
2010). Bolocat identifies sources based on their signifi¬ 
cance with respect to the local estimate of the noise in 
a map, subdividing regions with high significance into 
individual sources based on local maxima inside that re¬ 
gion. Each pixel is assigned to an individual source us¬ 
ing a seeded watershed, in a similar way to other source 
extraction algorithms ( Clumpfind , Williams et al. 1994; 
SExtractor , Bertin & Arnouts 1996). The source identifi¬ 
cation process is determined by three parameters: P amp , 
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the amplitude of the signal compared to a local esti¬ 
mation of the noise a (noise estimation as ck, 5); Phase , 
the base level of emission at which the identified region 
is expanded; and the deblending parameter Pdeb, used 
in decomposing regions of emission with multiple local 
maxima. First, Bolocat masks all data above P a mp of 
a particular image, and extends a region to include all 
connected regions of emission above Pb ase = S) since 
areas with marginal significance adjacent to regions of 
emission are probably real. Second, each identified re¬ 
gion is examined and subdivided according to the level 
of contrast between local maxima. For each pair of local 
maxima in a region, the amplitudes of emission, I\ and 
I 2 , are compared with the highest contour of emission 
containing the local maxima pair, I crit , and if any Ii or 
I 2 is less than Pdeb above I cr u , that local maximum is 
discarded as a subregion. Bolocat also avoids recover¬ 
ing sources (and sub-sources from the deblending pro¬ 
cess) with sizes less than the beam FWHM. A detailed 
description of the algorithm is presented in Rosolowsky 
et al. (2010). BGPS catalogs VI.0 and V2.1 used as 
parameters for the extraction of sources P arnp =2cr(l, 6), 
Phase—1 <t(/,6 ) and Pdeb— 0.5<r(Z, 6), with cr(l,b) the local 
noise estimation, as a function of Galactic coordinates, 
on the BGPS maps. 

We tested Bolocat in the SHARC-II maps varying 
the P amp , Phase and Pdeb parameters, checking which 
of them reproduce a “by eye” extraction of sources in 
those maps. Using the same parameters as BGPS, spu¬ 
rious low brightness sources were recovered across the 
maps, with several of them toward the noisy edges of the 
maps, and bright extended sources are sub-divided until 
the routine reaches small areas, comparable to the beam 
size of the maps. Better results were obtained consid¬ 
ering P arnp =3cr(a,5 ), P ba8e =la(a,S) and P de fr=lcr(<a, 5), 
and these values were used for all maps (inclusive those 
outside the BGPS coverage) and the following catalog of 
BPGS substructures at 350 /am. 

Figure 4 shows an example of one SHARC-II map 
and the corresponding image at 1.1 mm obtained from 
BGPS. The angular scales of SHARC-II and BGPS im¬ 
ages are 1.6 arc-seconds per pixel and and 7.2 arc-seconds 
per pixel, respectively. The better angular resolution of 
SHARC-II reveals substructures within the clumps iden¬ 
tified by the BGPS catalogue. Crowded fields and ex¬ 
tended clumps identified at 1.1 mm show small compo¬ 
nents and filaments when they are observed at 8.5" res¬ 
olution. 

4.3. Source recovery experiments 

We performed similar tests to the ones done for BGPS 
to assess catalog properties compared with the real dis¬ 
tribution of emission on the sky. For these tests, we in¬ 
serted sources in selected maps that do not have detected 
sources, or that have just a couple of compact identified 
sources. The selected maps are those of L030.15+0.00 (2 
identified sources), L078.92—0.19 (0), L079.62+0.49 (0), 
L079.ll—0.35 (1), L080.86+0.38 (2), Llll.26-0.77 (2), 
Llll.79+0.71 (0), L136.52+1.24 (2), and L137.69+1.46 
(0). The median rms noise of these maps is 667 mJy 
beam -1 , larger than 83% of the estimated noise in the 
whole set of maps. This set of maps will be consid¬ 
ered as a conservative representation of an emission free 
SHARC-II map. For those maps that already have iden¬ 


tified sources, we flagged the detection of those sources 
and any input sources near them, and therefore the re¬ 
covered properties of input sources will not be affected by 
the emission prior to the test. After the artificial Gaus¬ 
sian objects were added, the maps were reduced in the 
same standard way as the rest of the SHARC-II maps. 
Although Gaussian brightness profiles do not quite rep¬ 
resent all structures observed in the 350 /am maps, they 
are still good models of compact sources. Thus, the bulk 
of sources will be well represented in these tests. 

We examined the degree of completeness of the cata¬ 
log of substructures at different flux density limits. For 
the nine test maps, we input sources with FWHM of 9" 
and amplitudes uniformly ranging between 0.1cr rms and 
15cr rms , with cr rms the noise of each map. For a rrns = 667 
mJy beam -1 , this range correspond to sources with flux 
densities between 0.19 Jy to 28.5 Jy. Taking into account 
that the edges of the maps have a larger noise than the 
central area of the map with uniform noise, we consider 
the detection fraction of sources across the whole map, 
and sources inside a central area of 7' in diameter. Fig¬ 
ure 5 shows the results of the test of completeness. The 
vertical lines represent 1 to 6 times the average a rms 
value of the complete set of SHARC-II maps. For sources 
inside the central area of the map, the detection is >99% 
complete at the 6a rms level. For sources with amplitude 
less than 3cr rms , the detection fraction is less than 15%. 
For the rest of our test, we only considered sources inside 
a central area of radius 3.5' to avoid noisy edge effects. 

The next test we performed was examining the proper¬ 
ties recovered by the extracting algorithm for a distribu¬ 
tion of input sources. We compared input and extracted 
flux densities for objects with FWHM equal to 23" (2.7 
times the beam size of the 350 /am maps). Initially, we 
tested the flux recovery of artificial sources on apertures 
of 20" and 40", and the total integrated flux, in units 
of /iV. The integrated flux is estimated from the area 
of significant emission defined by the watershed extrac¬ 
tion. At an aperture smaller than the source size, we 
only recovered part of the total intensity as expected, 
and the flux recovery at a larger aperture as 40" is al¬ 
most equivalent to the integrated flux. All the sources 
extracted in our catalog have a recovered flux density 
less than 3000 /jlV (~ 300 Jy). Figure 6 shows results 
of this flux recovery test. For the integrated emission, 
calibration was done using the C beam conversion factor, 
and flux recovery at 20" and 40" was calibrated using 
C20 and C40, respectively. The image indicates that the 
40" aperture flux density could be underestimating the 
amount of flux recovered for input sources, with a differ¬ 
ence with respect to the curve of the integrated flux of 
~ 10%. That difference between integrated and aperture 
flux density recovery gets bigger for input sources with 
larger size. For example, for input sources with FWHM 
equal to 4 times the beam size (~30"), the difference be¬ 
tween a 40" aperture and the integrated flux is around 
35%. 

The size recovered with Bolocat for input sources with 
a peak signal of 50cr rms is shown in Figure 7. The sizes 
of small sources are well recovered by the algorithm, 
but they become underestimated for sources larger than 
~ 60". These results are similar to those found in the 
BGPS maps, where the radii of recovered sources become 
underestimated for radii > 200". The largest source that 
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we found in the 350 /a m maps has a size of ^48", with a 
major-axis of 65". 

Finally, we inspected how the algorithm decomposes 
sources into individual substructures. We tested how 
blended sources are identified as a single or as a couple 
of individual objects, considering pairs of input objects 
in the set of test maps. Each source of the pair of input 
fake sources have sizes between 14" and 34", and peak 
flux density of 50cr rms . Figure 8 shows the fraction of 
blended sources as a function of the separation between 
pairs of input sources. In general, sources at a distance 
less than the beam FWHM cannot be recovered individ¬ 
ually and they are assigned to the same source in the 
recovered catalog, and more than 50% of the pairs of in¬ 
put sources are recovered individually at distances larger 
than 20" (2.4 times the beam size). Figure 8 also shows 
that pairs of input sources with sizes less than 24" are 
better resolved at shorter distance than pairs with larger 
sizes than that. Pairs of sources are resolved more than 
50% at separations of 16" for small input sources, and at 
25" for larger sources. 

4.4. Catalog of sources in the 350 /am maps 

We have recovered from the Bolocat extraction 1386 
sources in the 350 /am maps. We name our sources in 
a similar way as described in the BGPS catalog, us¬ 
ing the peak position of the source in galactic coordi¬ 
nates: SHARCLGZZZ.ZZZZ db bb.bbbb. An additional digit 
with respect to the BGPS catalog was required to ac¬ 
count for the better resolution and smaller pixel size in 
the SHARC-II maps. Table 3 presents the properties re¬ 
covered for the extracted sources. The positions of the 
recovered sources in the SHARC-II maps are shown in 
Figure 25 in the Appendix. 

Figure 9 shows the distribution of recovered flux for 
these sources using three methods from our catalog: pho¬ 
tometry in 20" and 40" apertures, and the integrated 
flux. For the 20" aperture, the average and median value 
of the flux density are 23.15T1.59 Jy and 8.95 Jy, respec¬ 
tively. For the 40" aperture, the average and median 
values are 45.11T2.82 Jy and 16.63 Jy, resp. For the in¬ 
tegrated flux density, the average and median values are 
59.59T5.19 Jy and 12.06 Jy, resp. 

Figure 10 shows the deconvolved radii and the aspect 
ratio distribution. The deconvolved radius distribution 
extends to 46", and the catalog has 257 sources 19% 
of the total) with radii not resolved. These small sources 
are usually faint, with an average value of their peak 
signal-to-noise ratio of 4.4. The average and median of 
the radius distribution are 15.0" and 14.0", respectively. 
The aspect ratio is defined as a ma j /(Jmim and its average 
and median values are 1.53 and 1.45, respectively. A 
similar median aspect ratio is found in the BGPS V2.1 
catalog, and this could be a common feature of millimeter 
and submillimeter sources, or an artifact of the extraction 
algorithm. 

4.5. Correlation with BGPS sources 

We performed a spatial position matching analysis to 
get a correlation between sources found in the 350 /am 
maps and the BGPS V2.1 sources contained in them. 
Therefore, we excluded from this analysis the 3 maps of 
outer Galaxy regions outside the coverage area of BGPS. 


A simple criterion to match sources from catalogs at dif¬ 
ferent resolutions is to take the position of maximum in¬ 
tensity for sources of one catalog, and find which sources 
in the second catalog are 1 beam distant from that po¬ 
sition. For example, Contreras et al. (2013) performed 
this kind of matching between sources of the ATLAS- 
GAL survey (870 /im, 19.2" FWHM beam size) and 
BGPS VI.0 sources, considering that sources between 
these catalogs are associated if their peak positions are 
less than 40" in angular distance, finding around ^3000 
likely matches. One of the issues with this approach 
is that it does not yield a one-to-one association. This 
has special importance in catalogs based on sources ex¬ 
tracted from maps at different resolutions. Even more, 
Ginsburg et al. (2013) found that for the same 1.1 mm 
data at 33" FWHM beam resolution, re-processed V2.1 
maps with better spatial filtering recover in some cases 
more than the one source that previous VI.0 maps found, 
and in others V2.1 maps recover a single structure where 
VI.0 recovered many. Then, sources required more than 
simple matching at beam distance to compare different 
versions of the catalog. 

Considering the different beam size of the 350 /am maps 
(8.5") and 1.1 mm maps (33"), we expect to resolve 
extended millimeter sources in some cases into multi¬ 
ple smaller components, and we put the constraint that 
compact sources identified in the 350 /am maps must be 
associated with a single 1.1 mm source. For the follow¬ 
ing analysis, we will consider as a “parental clump” a 
1.1 mm source from the BGPS catalog harboring one or 
more associated matched structures on the 350 /am maps, 
which we will refer to as “high-resolution sources” or just 
simply as “substructures”. 

We made use of one of the sub-products of the Bolo¬ 
cat algorithm, the mask produced for each map in the 
identification of individual significant emission. These 
labelmask maps give information on which particular po¬ 
sition in a map contains emission, and associates that 
position with a single source from the catalog produced 
for that map. Our approach was taking the position of 
maximum intensity from the 350 /am sources, and if that 
peak is within the labelmask of a BGPS V2.1 source, it 
is considered associated. 

For the sample of 104 SHARC-II maps containing 
BGPS sources, there are 1374 high-resolution sources. 
Only 24 (2%) of those sources have their peak position 
not associated to a BGPS V2.1 source from the labelmask 
maps. The association between 350 /am high resolu¬ 
tion sources and 1.1 mm sources suffers from intrinsic, 
noise driven uncertainty in the peak flux density position, 
and the algorithm uncertainties in the labelmask area as¬ 
signed to a specific source. We improved the matching 
of sources considering also the labelmask maps of the 
350 /am emission, overlapping both masks and estimat¬ 
ing how much of the substructure mask area is associ¬ 
ated with the parental source in the BGPS labelmask 
area. Two of the 24 sources have most of their mask 
area associated with a parental source and therefore are 
re-considered as matched sources. For the sample 350 
/am sources with a match, there are 17 sources with their 
mask area shared between two parent sources, and for 
them we made a visual inspection to see the most likely 
correct parental association. Figure 11 shows a couple 
of examples for association between 350 /am and 1.1 mm 
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labelmask regions. 

We obtained a total of 1352 high-resolution sources 
associated with 349 different parental sources. The 
results of the previously described spatial match¬ 
ing of 350 fi m SHARC-II and 1.1 mm BGPS 
sources are shown in Figure 12. Half of the 
parental clumps have only one or two associated high- 
resolution sources, but toward several of the parental 
clumps we found a large multiplicity of sources, in¬ 
dicating possible fragmentation. The clumps with 
the largest multiplicity are BGPSv2_G213.705-12.603 
(34 substructures), BGPSv2_G034.256+00.154 (27), 
BGPSv2_G133.716+01.220 (25), BGPSv2_G000.014- 

00.017 (23), BGPSv2_G029.916-00.045 (22), BG- 

PSv2_G081.477+00.020 (20), and BGPSv2_G029.958- 
00.017 (19). There are 22 parental millimeter sources 
with more than 10 substructures, containing a total of 
350 high-resolution sources (26% of the total). 

These parental sources show some clear strong compact 
sources in them, but also some elongated, filamentary 
structures are present. We explored whether the large 
number of substructures in parental sources is due to 
spurious low brightness sources, and/or algorithm frag¬ 
mentation of large sources. We consider the peak signal- 
to-noise ratio identified for each source by Bolocat , and 
compare the number of identified sources under differ¬ 
ent ranges of this ratio. Table 4 presents the number 
of substructures with amplitude P a mp (described in sec¬ 
tion 4.2) above 6a rms , 10a rms and 20a rms for the 22 
BGPS sources with large multiplicity. For these, 66% 
of 350 fi m sources have P a mp of 6a rms , 42% of 10a rms , 
and 19% of 20a rms . Figure 11 shows two examples of 
the distribution at different peak signal-to-noise levels 
for two of these parental sources. While compact, strong 
sources have P a mp above 10a rms , detections below this 
limit recover not only isolated low emission objects, but 
also filament-like features in chains of sources, and weak 
detections produced in the decomposition process of the 
algorithm. We conclude then that sources with P a mp 
above 10a rms (“compact substructure”) are more related 
with possible dense core-like structures in the interior of 
millimeter clumps, and sources below this limit (“faint 
substructure”) trace a more diffuse medium in addition 
to weak sources. For the total number of high-resolution 
sources recovered from the SHARC-II maps, only 437 
(~32%) are considered compact substructures. Column 
14 in Table 3 indicates if a source recovered in the 350 
/i m maps is “compact” or “faint”. 

4.6. Flux densities of BGPS sources at 350 /am 

We determined the flux densities at 350 /am of BGPS 
V2.1 sources by measuring the flux density in an aper¬ 
ture of 40" on our sample of SHARC-II maps. To match 
the 33" FWHM effective BGPS beam at 1.1 mm, the 
350 fi m images, (1.618" per pixel) were convolved with 
gaussian kernels (using IRAF task GAUSS) having a 
= 13.5" (8.33 pixels), i.e., convolving with a gaussian 
FWHM= ay/sin 2 = 31.8". 

From the total of 619 BGPS V2.1 sources in the 350 
/am maps, there are 82 sources that appear in more than 
one map; we preferred the flux determination from the 
maps where those sources suffer fewer sampling artifacts, 
such as high noise due to map edge proximity, or neg¬ 


ative bowls around the source. In case there are not 
evident problems in the source flux estimation, we just 
considered the averaged result of the flux density. In 
addition, there are 45 BGPS sources with negative inte¬ 
grated fluxes on the 350 /am maps due to proximity to 
noisy edges or negative bowls surrounding areas of strong 
emission, and therefore these sources are not considered 
as reliable measurements. 

Results of the estimated values of 350 /am continuum 
emission from BGPS sources are presented in Table 5. 
Column 1 gives the name of the source in the V2.1 cat¬ 
alog. Column 2 gives the flux density integrated in an 
aperture of 40" centered on the peak position of the 1.1 
mm source. The values of the flux density are corrected 
by the factor 1.46 suggested by Aguirre et al. (2011) for 
the 40" aperture flux obtained in BGPS catalog. Column 
3 gives the integrated flux density for aperture photom¬ 
etry in the same previous position, but this time in the 
350 /am convolved maps. 

4.7. Comparison between SHARC-II and Herschel 
images 

From Molinari et al. (2010), we used the Hi-GAL 
2° x 2° images at 350 /am obtained during science demon¬ 
stration phase centered toward l = 30°, 5 = 0° and 
compared the results of flux recovery between the Her¬ 
schel /SPIRE image and the 10 SHARC-II maps con¬ 
tained in that area. The Herschel image has a FWHM 
beam resolution of 24.9", and details in the reduction 
process are found in detail in Traficante et al. (2011). 
Figure 26 in the Appendix shows 1.1 mm BGPS images, 
and their corresponding Hi-GAL images at 350 /am to¬ 
ward l = 30°, 5 = 0°, for regions mapped with SHARC- 
II. While the 1.1 mm images in this figure show the po¬ 
sition of the sources from the BGPS V2.1 catalog, the 
350 /am images show in blue (< 10<r) or red (> 10a) the 
objects identified in our high-resolution source catalog. 
There are 213 of these substructures found, with only 
one left out of the following analysis due to noisy edge 
effects. 

First, we compared the flux obtained on the sources 
recovered directly by Bolocat in this set of SHARC-II im¬ 
ages, with emission measured in the Herschel map. Up¬ 
per panel of Figure 13 shows the ratio between the den¬ 
sity flux at 40" aperture Fsharc and the recovered flux 
at the same aperture, centered on the peak position of 
the high-resolution sources, obtained from the Herschel 
image, FHerschei • For faint substructures, the flux ratio is 
in general below 0.3, with a cut around FHerschei ~ 100 
Jy, which suggests that the Herschel image is recover¬ 
ing additional flux from diffuse, large scale emission that 
SHARC-II maps do not recover. 

Uniform background emission ranging between 2700 
- 3300 MJy sr -1 in the Herschel image would give a 
flux density into a 40" aperture of 80 - 100 Jy. Fig¬ 
ure 14 shows one of the SHARC-II maps toward l = 30°, 
L029.95-0.05, and the Herschel emission in the same re¬ 
gion. This figure also presents the SHARC-II map con¬ 
volved to a resolution of 24.9", to match the beam size of 
Hi-GAL data at the same wavelength. The emission level 
at 17a in the Herschel data corresponds to ^2790 MJy 
sr -1 , and the morphology of emission recovered above 
that contour level in the Herschel map is similar to the 
structures detected above 3a level in the SHARC-II im- 
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ages. A “background emission” level for the Herschel 
maps would lie at about 2460 MJy sr _1 , corresponding 
to the 15cr level (white contour in Fig. 14), and a 40" 
aperture integrated flux density of ^73 Jy. The beam- 
matched SHARC-II image shows that the la emission 
level has a good resemblance with that background emis¬ 
sion level. 

Inspecting Figure 26, most of the faint substructures 
do not look like well-defined entities at a 10 cr contour, 
but instead they seem to be immersed in background 
extended emission. For the compact substructures re¬ 
covered in the SHARC-II maps, the average and me¬ 
dian values for the flux ratio are 0.56±0.03 and 0.53, 
respectively. For those high-emission compact sources 
with peak signal-to-noise above 50, represented by green 
points in Figure 13, flux recovery in the SHARC-II maps 
and the Herschel map are nearly the same (average flux 
ratio of 1 . 00 ± 0 . 11 ). 

We performed a similar analysis for the 102 BGPS 
sources found in this region. We first convolved the 350 
/a m Herschel image to match the 33" resolution of Bolo- 
cam, and compare later the 40" aperture flux emission 
obtained toward the peak position of BGPS sources. Re¬ 
sults are shown in bottom panel of Figure 13. Errors 
in the flux ratio consider conservative uncertainties of 
20% in the flux from the Herschel image, and 30% in 
fluxes from SHARC-II data. Sources without any sub¬ 
structure associated have low mm emission and flux ratio 
FsHARc/F H erschei on average lower than 0.1. For those 
1.1 mm parental clumps associated with compact sub¬ 
structures (shown in red in the figure), the average and 
median values of the flux ratio are 0.76±0.08 and 0.69, 
respectively. 

We emphasize that the differences arise from the bet¬ 
ter angular resolution of SHARC-II compared to Her¬ 
schel /SPIRE data, and not properly to the background 
subtraction; SHARC-II images are able to better disen¬ 
tangle the very dense and inner parts of clumps from the 
surrounding diffuse emission. 

5. ANALYSIS 

5.1. Temperature determination 

The data at 350 /a m and 1.1 mm can be used to define a 
color temperature. We will refer to this quantity through 
this analysis as the “dust temperature”, although vari¬ 
ations in temperature across and along the line of sight 
make this color temperature only a rough guide to the 
actual dust temperatures. 

The observed intensity of continuum emission is given 
by S v = f2&eam[l - e~ T ^ d ] B v (T d ), where Tt hea m is the 
beam solid angle (33" beam) for both the 1.1 mm data 
and the convolved, beam matched, 350 /am observations, 
r u d is the optical depth of the emitting dust at each fre¬ 
quency and B v (Td) is the Planck function at the dust 
temperature XJ. The optical depth in the submm to mm 
regime is considered proportional to . In the optically 
thin limit, the temperature can be estimated according 
to: 

D _ SsBO^m _ mm/kT d ) - 1] 

J-l' C* 3 I Q • y -L J 

bl.lmm Z'l. 1 mm[exp(/lZ'350Mm/&T , d ) - 1] 

£ 350^777 and aSi. 1777777 are the integrated flux densities ob¬ 


tained in 40" apertures in the convolved SHARC-II maps 
and in the BGPS maps, respectively. The distribution of 
the ratio £ 350^777 /£i.i mm for the 574 V2.1 sources with 
reliable flux values is shown in Figure 15. The ratio dis¬ 
tribution has a maximum of 181.0, with a median value 
of 19.7. 

Equation 1 cannot be solved analytically for 7% The 
adjustment of the pair of parameters X^ and /3 has been 
a recurrent issue in the study of the interstellar medium 
and star formation. Frequently, a value of /3^2 in the 
millimeter to submillimeter regime is assumed, referring 
to studies of dust grains composed by graphite and sil¬ 
icate (Draine & Lee, 1984), but models with variations 
of ice mantles due to coagulation give lower values (e.g., 
Ossenkopf & Henning, 1994, estimated (3 = 1.8 for their 
OH5 model). Values near /3 ~1 may be appropriate 
for circumstellar disks (e.g., Beckwith & Sargent, 1991). 
Shirley et al. (2005) summarize different opacity models 
in the submillimeter-millimeter regime for low-mass pre- 
protostellar cores, with values of /3 between 1.3 and 2.3. 
For the present analysis, we obtained dust temperatures 
considering three different values of [3: 1.0, 1.7, and 2.0. 
Nevertheless, considering that most BGPS sources prop¬ 
erties are related with dust clumps of dense material, 
(3 = 1.7 should be the closest model of the true nature 
of these structures. Figure 16 shows the fitted values for 
temperatures as a function of the ratio £ 350^777 / £ 1.1777777 
for the three models previously described. For /3 = 1.7, 
uncertainties in the determination of temperatures are 
shown in the figure, considering an error of the flux ratio 
of 10 %. 

The large uncertainty for high temperatures is un¬ 
avoidable with these data. At high temperatures (> 50 
K, see below), both 1.1 mm and 350 /am flux densities 
approach the Rayleigh-Jeans limit, and therefore the ra¬ 
tio becomes nearly constant, increasing significantly the 
uncertainty in the determination of X^. 

The computed temperatures using three different val¬ 
ues of (3 are presented in Table 5. Columns 4 and 5 gives 
the source color temperature determined from equation 
1 , using a emissivity index (3 = 1 . 0 , and the upper and 
lower limit for that temperature. Values of temperatures 
and upper and lower limits for a spectral index of (3 = 1.7 
and ^ = 2.0 are given in Columns 6-7 and 8-9, respec¬ 
tively. 

From equation 1, sources with flux ratio £ 350^777 / 
£ 1.1 777777 below 70.2 have temperatures lower than 1000 
K for /3 = 1.7. Only 32 V2.1 sources (~ 6% of sample) 
have estimated temperatures above 1000 K, most of them 
(24) found toward the Galactic Center region. In only 
one source, the large flux ratio is the consequence of low 
emission at 1.1 mm and comparable noise level, and al¬ 
most half of the rest have Si.i mm > 1 Jy, which included 
well-known regions such as G034.256±00.154, G029.958- 
00.017, G030.702-00.067, and the “Brick” IRDC (Long- 
more et al., 2012). These very high temperatures reflect 
flux density ratios near those expected in the Rayleigh- 
Jeans limit, and the temperatures are not constrained by 
our data. A centrally heated region will produce strong 
350 /am emission; using a ratio with the 1.1 mm emission, 
which traces a larger region of cold dust, can produce an 
artificially high color temperature. 

At a limiting temperature of 50 K, the corresponding 
flux ratio is 52.7 for /3 = 1.7, and a variation of 10% of 
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this ratio will give an increase of 40% in the estimated 
temperature. Therefore, to avoid large uncertainties we 
consider as a “good fit” a determined temperature less 
or equal than 50 K, and a similar limiting temperature 
value is found for (3 = 1.0 and (3 = 2.0. There are 512 
sources in our sample below this limit for (3 = 1.7, and 
30 sources with temperatures between 50 and 290 K with 
quite large uncertainties. 

Figure 17 shows the temperature distribution for the 
512 sources, assuming (3 = 1.7, with good fits for tem¬ 
perature. The median value is 13.3 K, with an average 
of 16.3±0.4 K. The difference between the median and 
the average is explained by the positive skewness of the 
distribution. The figure also shows the distribution of 
sources with different models for the spectral index. Con¬ 
sidering a value of (3 = 1.0, the sample of V2.1 sources 
tends to have larger values of temperature, and then the 
number of these sources with T d < 50 K is 322, with a 
median value of 17.3 K and an average of 19.6±0.5 K. 
For [3 = 2.0, the distribution of fitted temperatures has 
lower values, with 543 sources below 50 K, and with me¬ 
dian and average temperatures of 11.2 K and 13.3±0.3 
K, respectively. 

Figure 18 plots temperature versus flux density 
at 1.1 mm for a fixed value (3 = 1.7. For faint 
mm sources, uncertainties in the flux dominate the 
uncertainties in the temperature. In contrast, tem¬ 
peratures above 30-40 K are uncertain because of 
the weak constraints as the Rayleigh-Jeans limit is 
approached. The figure also shows the presence of 
8 sources with high flux density (> 4 Jy at 1 mm) 
and low estimated temperatures (< 20 K). These are 
BGPSv2_G359.946-00.045 on the Galactic center region, 
BGPSv2_G075.834±00.400 and BGPSv2_G076.384- 
00.622 on Cygnus-X, BGPSv2_G133.716±01.220 on W3 
main, BGPSv2_G031.411±00.307, and the sources BG- 
PSv2_G206.557-16.361 and BGPSv2_G206.534-16.356 
on NGC 2024, toward the Orion B South molecular 
cloud. In the case of BGPSv2_G133.716±01.220, 
corresponding to W3 East, Rivera-Ingraham et al. 
(2013) recovered a higher temperature 30 K for 
f3 = 2.0) for this source using Herschel data from 
the HOBYS program (Motte et al., 2010), difference 
that is explained in part for the underestimation of 
the flux at 350 fi m due to a shift on the peak of the 
BGPS source (and therefore the integration aperture), 
and the peak of the source on the SHARC-II map. In 
fact, there is a better agreement on the sources W3 
West (BGPSv2_G133.698±01.216; T d = 28 K, /3 = 2.0) 
and W3 SE (BGPSv2_G133.747±01.198; T d = 25 K, 
/3 = 2.0), and the values estimated by Rivera-Ingraham 
et al. for these clumps (^25 K and ^21 K, respectively). 

In the case of the clump associated to the G31.41±0.31 
hot core (Cesaroni et al., 1994), an isothermal tempera¬ 
ture of T = 29 K was estimated for this source by Mueller 
et al. (2002) from continuum data modeling, which could 
indicate that the low temperature is the result of unreli¬ 
ability on the measurement of derived fluxes. 

For the two clumps on the Orion B cloud, Johnstone et 
al. (2006) considered for both sources a temperature of 50 
K, though they were unable to obtain color temperatures 
for them from SCUBA 450 /im and 850 /am continuum 
data due to high flux ratios, inconsistent with modified 
blackbody emission. 


There are other 22 sources with Si.imm > 1 Jy and 
T d < 15 K. These sources could be dense, prestellar 
clumps, but further observations and analysis, particu¬ 
larly in molecular line observations of high-density trac¬ 
ers, are necessary to test their nature. 

5.2. Comparison of temperatures with other surveys 

We found in general a good agreement between our 
determination of temperatures and estimations of tem¬ 
perature from other surveys of high-mass star forming 
regions. In this section, we describe details of these com¬ 
parisons. 

Faundez et al. (2004) performed an analysis of physical 
properties of 146 continuum emission structures detected 
at 1.2 mm with SIMBA toward IRAS sources associated 
with CS(2—1) detections. The temperature and spectral 
index of these sources were obtained from SED fitting 
using additional IRAS four band fluxes, obtaining values 
for (3 ranging between 1.5 and 2.5, with average tem¬ 
perature for their sample of ^32 K. Temperatures found 
for their millimeter structures are larger than ours be¬ 
cause their sample was biased toward IRAS sources that 
are intrinsically warm. Their sample has strong emis¬ 
sion at 1.2 mm, with all their sources (except one) above 
1 Jy and with an average flux density of 16.3±3.2 Jy. 
Considering the emission at 1.1 mm from our sources 
that are, within uncertainties and corrections, compa¬ 
rable with their sample, the average temperature for the 
124 V2.1 sources with 40" aperture flux at 1.1 mm above 
1 Jy is ^38 K for /3 = 1.7. Similarly, Mueller et al. 
(2002) modeled SEDs of a group of 51 massive star form¬ 
ing clumps, originally selected by having water masers, 
and found that a temperature of 29 ± 9 K best charac¬ 
terized the sources. Cross-matching their sample with 
BGPS sources inside SHARC-II maps, considering one 
beam (33") in distance associations between samples, we 
found 13 BGPS clumps with good determination of dust 
temperature associated with the sources in Mueller et 
al. (2002). These sources have average and median val¬ 
ues for their 40" aperture integrated flux of 8.5 Jy and 
7 Jy, respectively, with only one source with integrated 
flux below 1 Jy. One of these sources corresponds to 
BGPSv2_G031.411±00.307, for which its fitted dust tem¬ 
perature is considered underestimated (see Section 5.1). 
Excluding this source, the fitted dust temperatures of 
the rest of the clumps have an average of 27 K, and a 
median of 29 K. Thus, the parameters that we found in 
sources with strong (>1 Jy) integrated emission at 1.1 
mm are similar to those found with SED fitting, giving 
some confidence to our values based on only two wave¬ 
lengths. In addition, BGPS clumps associated with mas¬ 
sive star forming regions generally have strong emission 
at mm and submm wavelengths, with warm tempera¬ 
tures ^30 K, lying above 89% of the distribution of fitted 
temperatures of BGPS sources with T<50 K. 

Battersby et al. (2011) used Hi-GAL 70 to 500 /am 
data to identify significant dust continuum emission in 
two 2°x2° regions centered at l = 30° and l = 59°, and 
obtained temperature and column density maps for them 
by fitting spectral energy distributions for each pixel. 
While the temperature is a free parameter for their fitting 
model, the spectral index is fixed to 1.75 following Os- 
senkopf & Henning (1994). From our sample, we have 91 
BGPS V2.1 sources with determined temperatures in the 
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I = 30° map, and none in the l = 59° map. Consider¬ 
ing Tsharc the fitted temperature obtained from equa¬ 
tion 1 , the median and average temperature for those 
sources from our analysis is 18 K and 344=4 K for /3 = 1.7. 
We estimated a value of the temperature of these sources 
from the maps presented by Battersby et ah, TnerscheU 
averaging the temperature values of each pixel on an 
aperture of 40" centered on the peak position of each 
source. Figure 19 shows the ratio T S harg / T Hersche i 
as a function of the flux density at 1.1 mm for the sam¬ 
ple of BGPS sources. The average ratio from this figure, 
without considering those sources with Tsharc > 50 K 
for which the uncertainties are much larger, is 0.83T0.05, 
while the weighted mean is 0.50T0.01. 

To identify dense material associated with the forma¬ 
tion of massive stars and clusters, Battersby et al. (2011) 
performed iterative removal of cirrus cloud emission on 
Hi-GAL images. The spectral emission distribution fit¬ 
ted for each point across the map results in a median 
value for temperature of T=23 K for a spectral index of 
/? =1.7 on l = 30°. The region considered as source 
emission is roughly the contour at 6.5cr (^1070 MJy 
sr -1 ) in the 350 fi m SPIRE band, and hence SHARC- 

II maps still miss extended emission that is considered 
to be source emission by Battersby et al.. From their 
temperature map, we estimated a temperature for this 
extended background emission, masking the region above 
15a from the 350 /im SPIRE map, and also masking cir¬ 
cular areas of 80" diameter centered on each BGPS V 2.1 
source on the map. We found a temperature of 25 K, 
for a fixed /? =1.75 for the emission considered for SED 
fitting by Battersby et al. but filtered out in SHARC-II 
maps. 

We tested then the estimation of color temperature on 
compact sources, if additional background emission were 
considered. A simple model of a source emitting uni¬ 
formly as a modified blackbody, characterized by a tem¬ 
perature T source and a spectral index /3 source , combined 
with background emission, also modeled as a modified 
blackbody, with Tback and fiback- The measured fluxes 
of the composite emission, 6350 ^™ and Si.i mm , corre¬ 
spond to the sum of individual contributions at those 
wavelengths, and following equation 1 , differences be¬ 
tween Tsource £md Tback, £md (3source £md (3back-) will re¬ 
cover a new set of measured values Tm and Pm- For a 
field near l = 30°, where the very dense structures are 
likely spread over diffuse emission with Tb ac k = 25 K 
and pback = 1-7, we simulated the recovery of properties 
for input sources with a range of parameters character¬ 
istic of clump-like features, 5 < T SOU rce < 50 K, and 
1.5 (3source ^ 2.5. The ratio Tsource fTsource-\-bg of the 

input to the determined temperatures, as a function of 
the measured 1 . 1 mm flux, are shown in the lower panel of 
Fig 19. The results are qualitatively similar to the com¬ 
parison between Tsharc and THerschei, with low-flux, 
cold sources having low ratios of temperatures; warm, 
strong mm emission clumps have ratios of temperatures 
slightly above 1 . 0 ; and with sources with temperatures 
of 10-20 K (the bulk of BGPS clumps in l = 30°) with 
temperature ratios in the range 0 .4-1.0. 

A more direct comparison can be performed with com¬ 
pact sources extracted on Herschel images using the CU- 
TEX method, designed for detection of sources on in¬ 


tense and highly variable fore/background (Molinari et 
al., 2011). Using a similar method as presented by Elia 
et al. (2013), a catalog of extracted compact sources and 
their associated physical parameters, using the five bands 
of Hi-GAL, have been obtained for the inner Galaxy 
(—70° < l < 67°) (D. Elia et al., in preparation). These 
sources are extracted from background emission and 
therefore considered to represent the very dense clump- 
core features across molecular clouds. From a simple 
spatial match performed to the 102 BGPS sources to¬ 
ward l = 30° that have estimated color-temperatures, 
76 of these clumps have a counterpart Hi-GAL com¬ 
pact source, and preliminary results produce much better 
agreement with the temperatures derived in this paper 
from SHARC and BGPS data; the ratio between these 
temperatures have a mean of 1.05 and a median of 0.97. 
In both cases, a fixed value of the spectral index /? = 2.0 
was considered. Several considerations need to be as¬ 
sessed, such as different beam resolutions between obser¬ 
vation bands, and a more complete analysis between the 
presented sample of SHARC-II maps and properties of 
extracted sources with Hi-GAL, including substructures 
and association with emission at 70 /im, will be the sub¬ 
ject of a forthcoming work. 

5.3. Comparison between gas temperature from NH% 
and color temperature 

Figure 20 shows the ratio of gas to dust temperatures 
for a sample of BGPS sources. The gas temperatures 
were determined from NH 3 observations (Dunham et al., 
2010; Wienen et al., 2012 , and BGPS team in prep.). 
For those observations not centered on peak positions 
of BGPS sources, we associated the NH 3 detections to 
their closest BGPS source within one beam distance. 
Only the 49 clumps with dust temperatures below 50 K 
and gas temperatures below 30 K are plotted. The ratio 
T gas )Td U st has a weighted average of 0.88 and a median 
of 0.76. There are few sources with equal temperatures, 
but instead a group with lower T gas and a group with 
higher T gas . 

5.4. Comparison with tracers of dense gas 

Schlingman et al. (2011) presented spectroscopic ob¬ 
servations of HCO + and N 2 H + (3-2) for 1882 sources 
from BGPS. The selection of sources was made be¬ 
tween Galactic longitudes 10 ° < l < 90°, considering 
all sources with integrated flux in a 40" aperture greater 
than Si.i mm ~0.4 Jy, and from bins of logarithmically 
spaced flux of 0.1 for sources with Si.i mm = 0.1 — 0.4 
Jy. A new, complete version of this catalog by Shirley 
et al. (2013) includes all 6194 sources in the BGPS VI .0 
catalog between 7.5° < l < 194° of this pair of molec¬ 
ular transitions. After considering a one-to-one match 
between sources from VI.0 and V 2.1 catalogs, with a spa¬ 
tial matching of one beam of distance (33"), we checked 
which of those sources are present in our sample of 
SHARC-II maps. There are 300 V 2.1 sources with spec¬ 
troscopy data, 250 of them with HCO + detections and 
estimated dust temperatures from our analysis, and 199 
with N 2 H + detections with dust temperatures. Figure 21 
shows the peak main beam temperature versus our fitted 
temperature for that group of sources. The points scatter 
widely about the plot. Considering only the sources with 
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dust temperatures less than 50 K (169 for N 2 H + and 223 
for HCO + ), the Spearman’s rank correlation coefficient 

of the distributions of log (T^) and log (Tdust) are 0-51 for 
sources with N 2 H + detections, and 0.33 for sources with 
HCO + detections. We tested the null hypothesis of no 
significant correlation between log(X^) and log (Tdust) 
for these distributions. For a sample size n, the Student’s 
t value for the Pearson correlation coefficient r(n) is es¬ 
timated from t = r/s r , with s r = yf(\ — r 2 )/(n — 2). 
For a two-tailed probability value p = 0.05, our sam¬ 
ple requires t — 2, corresponding to r(169) = 0.151 and 
r( 223) = 0.131. The estimated values for the correlation 
coefficients of our sample are r = 0.508 and r = 0.332 for 
N 2 H+ and HCO + , respectively, rejecting the null hypoth¬ 
esis at a 95% confidence. Therefore, although the corre¬ 
lation is weak, in both cases it is statistically significant 
for the detected peak main beam temperature and the 
fitted dust temperature. In the case of N 2 H + , the corre¬ 
lation can be fitted as log(X^) = Bxlog(T^) + A, with 
fitted values A=—1.55 and B=1.04±0.13. For HCO + , 
the fitted parameters are A=—0.59 and B=0.52±0.08. 

Shirley et al. (2013) suggested that BGPS sources with 
a higher N 2 H + /HCO + ratio could be clumps with dense 
core-like structures in them, as a consequence of cold, 
CO-depleted gas within those cores. Nevertheless, we 
did not find a clear correlation between the ratio of high 
density tracers N 2 H + /HCO + , and number of substruc¬ 
tures or estimated dust temperature in our source sam¬ 
ple. Parental clumps with large multiplicity of substruc¬ 
tures and high N 2 H + /HCO + ratios represent good can¬ 
didates for future interferometric observations to study 
core properties in regions of clustered massive-star for¬ 
mation. 


6. DISCUSSION 

6.1. Mass and surface density of substructures 
We use the following equation to calculate the mass: 


M d 


n v B v (T d ) 


( 2 ) 


where S v is the flux density, D is the distance to the 
source, and k v is the dust mass opacity coefficient. 

We use a dust opacity of 1.14 cm 2 g -1 of dust at 1.1 
mm (Ossenkopf & Henning, 1994) for a model with thin 
ice mantles (OH5). Assuming a dust-to-gas mass ratio 
of R = M d /M g = 0.01, the mass of gas and dust will be 
given by: 


Afl.lmm = 14.26 



2 


1.14 cnv g 




(3) 


From our sample of 512 BGPS sources with temper¬ 
atures < 50 K, we can give some estimates of how the 
values of masses will change when a fixed temperature is 
used. The average ratio between the clump mass M(T d ), 
with T d the fitted dust temperature from continuum ob¬ 
servations between 350 pm and 1.1 mm, and the mass 
estimated at a fixed temperature, M(Tf ix ), is obtained 
from: 


/ M(T d ) \ = / exp(13.01/T d ) -1 \ 

\ M(T fix ) / \ exp(13.01/T /ix ) - 1 / 1 ’ 

Figure 22 shows the results of the above equation for 
different values of Tfix- It is common to use for simplic¬ 
ity a fixed value of 20 K in determination of masses from 
millimeter surveys of Galactic star forming regions (e.g. 
Motte et ah, 2007; Schlingman et ah, 2011; Urquhart et 
ah, 2013). Using Tf ix = 20 K, the average value of the 
ratio is 2.43T0.11, with a median value of 1.83. There¬ 
fore, assuming a fiducial value for the temperature of 20 
K for our sample of clumps, the estimation of masses 
would be underestimated on average, with respect to the 
mass obtained from our analysis of color temperatures. 

From Ellsworth-Bowers et al. (2013, 2015), kinematic 
distances were determined for a large number of BGPS 
sources. Only a fraction of our sample of BGPS clumps 
have estimated distances, so we cannot calculate masses 
for the entire group. However, we can estimate the con¬ 
tribution of masses from substructures detected at 350 
um with respect to the parental clump, because both 
substructures and clump are located at the same dis¬ 
tance. We use the flux density at 350 pm to compute 
masses of the substructures. Similar to equation 3, but 
with appropriate change in the exponential and with a 
dust opacity of 11 cm 2 g _1 (OH5), the mass measured at 
350 pm is as follows: 



From our sample of 349 BGPS parental clumps asso¬ 
ciated with 1352 high-resolution sources, the fraction of 
the total mass estimated for structures at 350 pm with 
respect to the mass of the parental 1.1 mm clump is de¬ 
termined by: 


T, (-^substructures) 


M ( 


clump 


3.75 x 10" 3 G 350 *"" ) 

V ^ 1.1mm ) 

( exp(41.14/Trf) — 1 \ 

X Ve x P( 13 - 01 / r d) - 1 J 


,( 6 ) 


where Sf lrnm is the emission of the parental BGPS 
clump and S^ 50jLim represent the total emission of sub¬ 
structures inside that clump. We consider the same dust 
temperature T dl fitted from equation 1, for substructures 
and their parental clump. Here the integrated emission 
recovered by Bolocat for each clump is taken into ac¬ 
count, which estimates the total flux for the area delim¬ 
ited by the labelmask region (see Section 4.5). In that 
way, the contribution of each substructure will be associ¬ 
ated with a single parental clump, with both considered 
at the same dust temperature. 

Figure 23 shows the ratio between masses obtained 
for high-resolution sources at 350 pm and their parental 
BGPS clump, as a function of the dust temperature 
and as a function of the integrated flux at 1.1 mm. 
Black points represent the mass ratio considering the 
flux contribution of all high-resolution sources, and red 
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points only consider emission from strong, compact 
sources. The average and median values for the mass 
ratio {^substructures ) / M c i urnp when all substructures 
are considered, are 0.48 and 0.40, respectively, with a 
mean weighted by errors of 0.22T0.01. Parental clumps 
that are more massive and cold have smaller uncertain¬ 
ties, so the mean mass ratio weighted inversely by un¬ 
certainty is smaller. When only compact substructures 
(see Section 4.5) are considered, the average and median 
of the distribution are 0.39 and 0.37, respectively, with 
a mean weighted by errors of 0.19T0.01. These results 
for the mass fraction are larger than the results from 
Ragan et al. (2013), who fit SEDs and obtained a value 
of ~ 14% for the mass contribution of core-like features 
observed at a resolution of 7.8" compared to the total 
mass of their harboring IRDCs. Possible explanations of 
this difference are our larger sample of sources (Ragan et 
al. consider 11 IRDCs, harboring a total of 83 substruc¬ 
tures), and also that, contrary to BGPS sources, IRDCs 
are not limited to compact structures. Besides, our sam¬ 
ple of clumps has a broader range of evolutionary stages, 
not necessarily limited to the coldest, initial stages of 
massive star formation, represented by the IRDC phase. 

Assuming a uniform temperature will also introduce er¬ 
rors in the mass determination. The parental molecular 
clump will be affected by both embedded stars and the 
external radiation field, and it is likely that substructures 
will have a different temperature from their surrounding 
medium. In addition, SHARC-II and Bolocam are sensi¬ 
tive to different spatial scale emission (see Section 4.3), 
and therefore the recovered flux at 1.1 mm could include 
faint, extended background emission not detected in 350 
/a m continuum emission. A better analysis of the mass 
ratio between parental clumps and their internal struc¬ 
tures (compact, core-type sources, and faint, filamentary 
emission) will require a comparison between SHARC-II 
and Herschel maps at 350 /a m, or modeling these sources 
using radiative transfer codes. 

Figure 24 shows the surface density E 3 5 o Mm estimated 
for the substructures found in the 350 /a m map, as a 
function of the surface density of their parental clump, 
We considered the deconvolved radius of re¬ 
solved substructures at 8.5"and resolved clumps at 33". 
Identified substructures have in general larger surface 
densities than their parental BGPS source. For the sam¬ 
ple of spatially resolved sources, there are 418 compact 
substructures, 286 of which have estimated temperatures 
< 50 K. These are represented by red points in the fig¬ 
ure. A linear fit of the compact sources indicate that they 
have £ 350^771 ~ 2.2x£i.i mm . Eighteen of them ( 6 %) have 
S 350 ^m > 1.0 g cm -2 , a theoretically predicted threshold 
for the formation of massive stars (Krumholz & McKee, 
2008). More recently a lower value of 0.2 g cm -2 was 
proposed by Butler & Tan (2012) as the initial condition 
of massive star formation in the Galaxy; 88 (31%) of our 
sources meet that criterion. 

Thus, for our sample of high-resolution structures ob¬ 
served at 350 /im, there are probably many dense cores, 
whose physical parameters suggest on-going or imminent 
formation of massive stars. 

6 . 2 . Uncertainties on determined physical parameters of 

clumps 


Following equation 1, the presence of additional inte¬ 
grated flux densities observed at 350 /am and 1.1 mm not 
related to the dust emission of a given clump-like source, 
such as emission from the surrounding environment of 
the compact source, or contamination from molecular 
lines on those bands, can introduce dissimilarities be¬ 
tween the recovered and real flux ratio and therefore 
in subsequent derived physical parameters. As seen in 
Section 5.2, the temperature estimated for a source em¬ 
bedded in a background with different temperature (or 
spectral index) may differ from the actual source temper¬ 
ature. While the data-reduction pipeline of the Hi-GAL 
survey is coded to remove low frequency noise, which im¬ 
ply that the 350 /am Herschel images recover faint diffuse 
emission (Traficante et al., 2011), SHARC-II maps filter 
out most of that emission and recover the very dense, 
compact structures, and BGPS maps are not sensitive 
to extended emission at scales larger than 2'-3' (Gins- 
burg et al., 2013). The precise amount of additional flux 
per band is determined by the temperature and spectral 
index of the background extended emission (see Fig. 19). 

In the case that observations of a particular clump re¬ 
cover some background emission in the 1.1 mm BGPS 
but not in the 350 /am SHARC-II maps, their overall ef¬ 
fect would consist of lower measured color temperature, 
and higher mass estimation. Contamination of CO(7 — 6 ) 
lines by the 350 /am emission is not expected for the bulk 
of clumps in this work. The J = 7 level of CO is 155 K 
above ground, and even though this line has been seen in 
sources with kinetic gas temperature of ^25 K (Krugel 
et al., 1989), emission is expected only in the warmer and 
denser clumps. No line contamination is possible at 1.1 
mm since the passband of Bolocam excludes the strong 
CO (2 - 1 ) line. 

To quantify these uncertainties, we consider two hypo¬ 
thetical cases of clumps with equal spectral index (3 = 1.7 
and temperatures T source of 20 K and 40 K, respectively. 
These sources have an associated value R SOU rce of the ra¬ 
tio of integrated flux densities observed at 350 /am and 
1.1 mm, and additional emission in these bands will pro¬ 
duce variation of the measured value of the flux ratio, 
Rm- For the case T source = 20 K, differences of less 
than 20 % between Rm and the intrinsic flux ratio of the 
source R SO urce, will imply a measured temperature Tm 
also 20% different from the source temperature T source . 
From equation 3, the uncertainty in the determination of 
the clump mass and surface density will be ^30%. 

For warmer sources, even small changes to the flux ra¬ 
tio will result on large uncertainties in the measured Tm- 
For T source = 40 K, differences of 10% between Rm and 
Rsource will result in differences of 30% of Tm with re¬ 
spect to T ] sourcei and clump mass and surface density will 
also have uncertainties ^30%. In the case of differences 
of 20 % between Rm and R SOU rce > measured temperatures 
could be 80% different from the real T source , and clump 
mass and surface density will have additional uncertain¬ 
ties of ^50%. 

7. SUMMARY AND CONCLUSIONS 

We presented a set of 107 SHARC-II 350 /am maps to¬ 
ward star forming clumps. Most of maps are observed 
toward sources from the BGPS V 2.1 catalog, with ad¬ 
ditional maps covering three outer Galaxy star forming 
regions. The maps have a resolution of 8.5", improving 
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upon the resolution of the BGPS maps at 1.1 mm (33 ") 
and revealing a population of substructures. 

We used Bolocat to identify significant emission in the 
SHARC-II maps, and we obtained a catalog of 1386 high- 
resolution sources, with average fluxes of 23.15±1.59 Jy 
and average radius of 15". We found that ^32% of these 
features have peak emission above 10 cr rms , and we con¬ 
sider these “compact substructures” as core-like sources 
immersed in parental clumps. Below this limit, recovered 
sources are called “faint substructures”, and they are re¬ 
lated with fragmentation of filamentary features, isolated 
low intensity sources, and residuals from the source ex¬ 
traction algorithm. 

We found 619 BGPS V 2.1 sources in our set of 350 fi m 
maps, only 56% of which are associated with 1352 high- 
resolution substructures. While most of the parental 
clumps have only one or two associated substructures, 
we found significant multiplicity in some of them, with 
22 clumps having more than 10 substructures in them. 

We compared the 350 //m continuum emission from 
SHARC-II maps with Herschel images from the Hi-GAL 
survey toward the l = 30° region at same wavelengths. 
Strong emission detected in SHARC-II maps is in gen¬ 
eral in good agreement with the Hershel fluxes, but faint 
sources from the SHARC-II maps are immersed in back¬ 
ground emission. 

We estimated temperatures for different models of the 
spectral index. For the 512 sources with “good fits” in 
the determination of temperatures (Td < 50 K), the me¬ 
dian and average values of their temperatures are 13.3 
K and 16.3±0.4 K, respectively, in reasonable agreement 
with other temperature determination methods based on 
SED analysis. 

Comparison to gas temperatures derived from NH 3 ob¬ 
servations reveals sources with ratios of Tg as /T(i U st both 
larger and smaller than unity, but with mean value of 
0 . 88 . 

Sources with N 2 H + and HCO + detections present 
weak correlations between their observed peak main 
beam temperature and their fitted dust tempera¬ 
ture. There is no clear correlation between the ratio 
T^(N 2 H+)/T^(HCO+) and the number of substruc¬ 
tures found on parental clumps. 

The fraction of the mass contained in substructures 
compared with the total parental clump mass has an 
average of 0.48, but this value decreases to 0.22 when 
weighted by uncertainties. Considering only the compact 
structures, likely to be core sources, the mass fraction be¬ 
tween substructures and parental clumps is 0.19, higher 
than other studies of high angular resolution at 350 fi m. 

Small, but significant fractions of substructures satisfy 
theoretical thresholds for massive star formation of 0.2 
(31%) or 1.0 ( 6 %) g cm -2 . 
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Table 1 

Average Calibration Factors 


Date 

Pfeeam 

(Jy beam -1 iV -1 ) 

C 20 

(Jy uV- 1 ) 

c 40 

(Jy uV- 1 ) 

P 33 " 

U 40 1 

(Jy uV- 1 ) 

2006 Jun 

9.84i0.68 

0.35i0.03 

0.28i0.02 

0.42i0.03 

2007 Jul 

7.01i0.21 

0.23i0.01 

0.20i0.01 

0.29i0.01 

2007 Oct 

8.08T0.32 

0.27i0.01 

0.23i0.01 

0.34i0.01 

2009 Sep 

6.05i0.25 

0.21i0.02 

0.17i0.01 

0.26i0.02 

2009 Dec 

5.89i0.11 

0.22i0.01 

0.17i0.01 

0.25i0.01 

2010 Jul a 

7.31i0.25 

0.24i0.01 

0.20i0.01 

0.30i0.01 

2010 Jul 24 

17.83i2.50 

0.67i0.06 

0.51i0.06 

0.77i0.09 

2010 Dec 

6.69i0.16 

0.23i0.01 

0.19i0.01 

0.28i0.01 

2011 Dec 

10.82T0.39 

0.38i0.01 

0.31i0.01 

0.46i0.02 

2012 Sep 

15.59±0.85 

0.52i0.03 

0.43i0.02 

0.65i0.04 


a Does not include calibration from July 24, 2010. 
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Table 2 

Observing Information 


Map 

Config. a 

Size 

Map Center 

R.A. (J2000) Dec. (J2000) 

Obs. Date 

1 a Noise b 
(mjy beam -1 ) 

L359.85+0.00 c 

A 

10 ' 

X 

10 ' 

17:45:15.356 

-29:03:48.64 

2006Jun 

994 

L000.00+0.00 c 

A 

10 ' 

X 

10 ' 

17:45:38.222 

-28:56:15.83 

2006Jun 

1453 

L000.15+0.00 c 

A 

10 ' 

X 

10 ' 

17:45:59.055 

-28:48:38.10 

2006Jun 

1438 

L000.30+0.00 c 

A 

10 ' 

X 

10 ' 

17:46:20.056 

-28:40:51.99 

2006Jun 

1086 

L029.95-0.05 

A 

10 ' 

X 

10 ' 

18:46:10.078 

-02:40:45.19 

2006Jun 

514 

L030.00+0.00 

A 

10 ' 

X 

10 ' 

18:46:05.764 

-02:36:40.25 

2006Jun 

462 

L030.15+0.00 

A 

10 ' 

X 

10 ' 

18:46:21.926 

-02:28:41.39 

2006Jun 

426 

L030.30+0.00 

A 

10 ' 

X 

10 ' 

18:46:38.373 

-02:20:39.29 

2006Jun 

496 

L030.45+0.00 

A 

10 ' 

X 

10 ' 

18:46:54.711 

-02:12:38.81 

2006Jun 

603 

L030.60+0.00 

A 

10 ' 

X 

10 ' 

18:47:11.158 

-02:04:35.09 

2006Jun 

888 

L030.70-0.07 

A 

10 ' 

X 

10 ' 

18:47:36.576 

-02:01:32.42 

2006Jun 

510 

L030.80-0.05 

A 

10 ' 

X 

10 ' 

18:47:43.784 

-01:55:09.17 

2006Jun 

529 

L030.88+0.13 

A 

10 ' 

X 

10 ' 

18:47:13.872 

-01:45:57.32 

2006Jun 

367 

L034.26+0.15 

A 

10 ' 

X 

10 ' 

18:53:19.759 

+01:14:36.57 

2006Jun 

395 

L076.16+0.09 

A 

10 ' 

X 

10 ' 

20:23:52.185 

+37:36:41.20 

2006Jun 

249 

L077.93+0.02 

A 

10 ' 

X 

10 ' 

20:29:23.675 

+39:00:58.88 

2006Jun 

389 

L078.14-0.32 

A 

10 ' 

X 

10 ' 

20:31:26.817 

+38:58:53.21 

2006Jun 

408 

L078.96+0.37 

A 

10 ' 

X 

10 ' 

20:31:03.924 

+40:02:54.24 

2006Jun 

445 

L079.28+0.30 

A 

10 ' 

X 

10 ' 

20:32:20.042 

+40:15:45.21 

2006Jun 

558 

L080.92-0.ll 

A 

10 ' 

X 

10 ' 

20:39:16.364 

+41:19:41.83 

2006Jun 

724 

L081.45+0.04 

A 

10 ' 

X 

10 ' 

20:40:23.200 

+41:50:33.25 

2006Jun 

274 

L081.68+0.54 

A 

10 ' 

X 

10 ' 

20:39:01.659 

+42:19:37.11 

2006Jun 

802 

L082.55+0.14 

A 

10 ' 

X 

10 ' 

20:43:35.922 

+42:46:07.20 

2006Jun 

255 

L075.76+0.41 

A 

10 ' 

X 

10 ' 

20:21:24.076 

+37:27:55.95 

2007Jul 

652 

L076.12-0.24 

B 

6 ' 

X 

6 ' 

20:25:07.547 

+37:23:14.34 

2007Jul 

388 

L076.35-0.58 

B 

6 ' 

X 

6 ' 

20:27:10.982 

+37:22:29.34 

2007Jul 

430 

L078.92-0.19 

A 

10 ' 

X 

10 ' 

20:33:18.191 

+39:41:24.34 

2007Jul 

689 

L078.17-0.31 

A 

10 ' 

X 

10 ' 

20:31:29.785 

+39:00:20.75 

2007Jul 

800 

L079.62+0.49 

A 

10 ' 

X 

10 ' 

20:32:36.698 

+40:39:16.34 

2007Jul 

1489 

L079.11-0.35 

A 

10 ' 

X 

10 ' 

20:34:33.491 

+39:44:46.95 

2007Jul 

1226 

L080.86+0.38 

A 

10 ' 

X 

10 ' 

20:37:00.649 

+41:34:35.34 

2007Jul 

1874 

L110.11+0.05 

C 

4' 

X 

4' 

23:05:11.371 

+60:14:41.67 

20070ct 

288 

Llll.62+0.38 

C 

4' 

X 

4' 

23:15:32.196 

+61:07:30.77 

20070ct 

342 

L134.28+0.86 

C 

4' 

X 

4' 

02:29:02.834 

+61:33:28.83 

20070ct 

275 

L134.83+1.31 

C 

4' 

X 

4' 

02:34:45.364 

+61:46:15.77 

20070ct 

296 

L136.38+2.27 

C 

4' 

X 

4' 

02:50:08.515 

+61:59:54.27 

20070ct 

257 

L136.83+1.07 

C 

4' 

X 

4' 

02:49:04.352 

+60:43:23.95 

20070ct 

250 

Llll.28-0.66 

A 

10 ' 

X 

10 ' 

23:16:04.681 

+60:02:06.00 

2009Sep-2009Dec 

286 

Llll.54+0.78 

A 

10 ' 

X 

10 ' 

23:13:44.302 

+61:28:10.18 

2009Sep 

496 

L133.71+1.21 

A 

10 ' 

X 

10 ' 

02:25:41.066 

+62:05:42.68 

2009Sep-2009Dec 

338 

L133.95+1.06 

A 

10 ' 

X 

10 ' 

02:27:03.912 

+61:52:14.05 

2009Sep-2009Dec 

338 

Llll.26-0.77 

A 

10 ' 

X 

10 ' 

23:16:11.144 

+59:55:27.74 

2009Dec 

538 

Llll.78+0.59 

A 

10 ' 

X 

10 ' 

23:16:13.500 

+61:22:51.01 

2009Dec 

412 

Llll.79+0.71 

D,A 

10 ' 

X 

10 ' 

23:15:52.220 

+61:30:02.26 

2009Dec-2012Sep 

307 

Llll.88+0.82 

D,A 

10 ' 

X 

10 ' 

23:16:15.419 

+61:37:42.69 

2009Dec 

473 

L136.52+1.24 

D 

10 ' 

X 

10 ' 

02:47:25.531 

+61:00:34.20 

2009Dec 

260 

L136.85+1.14 

E 

3' 

X 

3' 

02:49:28.155 

+60:47:02.96 

2009Dec 

337 

L136.95+1.09 

D 

10 ' 

X 

10 ' 

02:50:02.727 

+60:41:52.32 

2009Dec 

475 

L137.69+1.46 

H 

IB 

X 

IB 

02:56:47.527 

+60:41:21.90 

2009Dec 

667 

L138.30+1.56 

D,A 

10 ' 

X 

10 ' 

03:01:34.013 

+60:29:10.72 

2009Dec-2011Dec-2012Sep 

205 

L138.48+1.63 

D,A 

10 ' 

X 

10 ' 

03:03:09.052 

+60:27:39.61 

2009Dec-2011Dec-2012Sep 

266 

L173.14+2.36 

A 

10 ' 

X 

10 ' 

05:37:57.556 

+36:00:18.63 

2009Dec 

777 

L173.17+2.35 

G 

4' 

X 

4' 

05:37:59.420 

+35:58:27.54 

2009Dec 

320 

L173.47+2.43 

F 

2.5' 

X 

2.5' 

05:39:07.589 

+35:46:02.82 

2009Dec 

487 

L173.57+2.44 

F 

2.5' 

X 

2.5' 

05:39:24.823 

+35:40:55.71 

2009Dec 

421 

L173.62+2.81 

D 

10 ' 

X 

10 ' 

05:41:07.430 

+35:50:21.17 

2009Dec 

359 

L173.72+2.70 

G 

4' 

X 

4' 

05:40:52.684 

+35:41:45.31 

2009Dec 

387 

L173.76+2.67 

G 

4' 

X 

4' 

05:40:51.995 

+35:38:53.79 

2009Dec 

301 

L188.79+1.03 

A 

10 ' 

X 

10 ' 

06:09:06.309 

+21:50:45.79 

2009Dec 

259 

L188.95+0.88 

D,A 

10 ' 

X 

10 ' 

06:08:52.987 

+21:38:19.76 

2009Dec 

268 

L189.03+0.78 

D,A 

10 ' 

X 

10 ' 

06:08:39.856 

+21:31:11.38 

2009Dec 

290 

L189.12+0.64 

A 

10 ' 

X 

10 ' 

06:08:20.114 

+21:22:04.54 

2009Dec 

403 

L189.68+0.19 

D,A 

10 ' 

X 

10 ' 

06:07:47.653 

+20:39:28.06 

2009Dec 

155 

L189.85+0.39 

I 

18' 

X 

18' 

06:08:53.475 

+20:36:24.00 

2009Dec 

505 

L190.17+0.74 

A 

10 ' 

X 

10 ' 

06:10:51.620 

+20:29:49.63 

2009Dec 

184 

L192.60-0.16 

D 

10 ' 

X 

10 ' 

06:12:28.348 

+17:56:17.12 

2009Dec 

340 

L192.60-0.05 

D,A 

10 ' 

X 

10 ' 

06:12:53.279 

+17:59:28.55 

2009Dec-2010Dec 

271 

L192.72+0.04 

A 

10 ' 

X 

10 ' 

06:13:28.147 

+17:55:38.60 

2009Dec 

335 

L192.81+0.11 

D 

10 ' 

X 

10 ' 

06:13:55.135 

+17:53:09.02 

2009Dec 

328 

L192.98+0.14 

D 

10 ' 

X 

10 ' 

06:14:21.502 

+17:45:02.52 

2009Dec 

383 

L196.42-l.66 

A 

10 ' 

X 

10 ' 

06:14:36.606 

+13:52:04.40 

2009Dec 

470 

L203.23+2.06 

A 

10 ' 

X 

10 ' 

06:41:00.779 

+09:33:56.58 

2009Dec 

269 
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Map 

Config. a 

Size 

Map Center 

R.A. (J2000) Dec. (J2000) 

Obs. Date 

1 a Noise b 
(mjy beam -1 ) 

L203.35+2.03 

A 

10 ' 

X 

10 ' 

06:41:07.741 

+09:26:52.53 

2009Dec 

314 

L213.71-12.62 

D 

10 ' 

X 

10 ' 

06:07:42.663 

-06:23:27.15 

2009Dec 

421 

L217.37-0.07 

A 

10 ' 

X 

10 ' 

06:59:17.374 

-03:59:14.04 

2009Dec 

704 

L234.57+0.82 

D 

10 ' 

X 

10 ' 

07:35:28.482 

-18:45:34.36 

2009Dec 

1051 

L001.10-0.07 

A 

10 ' 

X 

10 ' 

17:48:28.655 

-28:01:44.77 

2010Jul 

849 

L023.31-0.26 

A 

10 ' 

X 

10 ' 

18:34:39.947 

-08:40:36.21 

2010Jul 

379 

L023.43-0.22 

A 

10 ' 

X 

10 ' 

18:34:45.046 

-08:32:55.27 

2010Jul 

375 

L024.50-0.08 

A 

10 ' 

X 

10 ' 

18:36:14.747 

-07:31:49.77 

2010Jul 

306 

L024.65-0.13 

A 

10 ' 

X 

10 ' 

18:36:44.447 

-07:25:23.71 

2010Jul 

357 

L024.78+0.12 

A 

10 ' 

X 

10 ' 

18:36:04.446 

-07:11:28.75 

2010Jul 

363 

L025.40-0.18 

A 

10 ' 

X 

10 ' 

18:38:16.246 

-06:47:09.20 

2010Jul 

482 

L030.61+0.16 

A 

10 ' 

X 

10 ' 

18:46:39.147 

-01:59:36.24 

2010Jul 

349 

L031.28+0.05 

A 

10 ' 

X 

10 ' 

18:48:14.651 

-01:26:51.80 

2010Jul 

303 

L081.11-0.16 

A 

10 ' 

X 

10 ' 

20:40:05.498 

+41:26:56.24 

2010Jul 

341 

L081.28+1.01 

A 

10 ' 

X 

10 ' 

20:35:39.600 

+42:17:38.79 

2010Jul 

250 

L081.39+0.73 

A 

10 ' 

X 

10 ' 

20:37:12.963 

+42:12:38.77 

2010Jul 

341 

L081.48+0.00 

A 

10 ' 

X 

10 ' 

20:40:37.896 

+41:50:26.77 

2010Jul 

353 

L081.76+0.60 

A 

10 ' 

X 

10 ' 

20:39:01.567 

+42:25:36.73 

2010Jul 

329 

L031.41+0.31 

A 

10 ' 

X 

10 ' 

18:47:34.152 

-01:12:46.30 

2010Jul24 

791 

L081.88+0.77 

A 

10 ' 

X 

10 ' 

20:38:39.559 

+42:37:37.68 

2010Jul24 

717 

L183.40-0.58 

A 

10 ' 

X 

10 ' 

05:51:16.720 

+25:43:32.50 

2010Dec 

255 

L189.78+0.33 

A 

10 ' 

X 

10 ' 

06:08:32.990 

+20:38:33.34 

2010Dec 

294 

L202.58+2.42 

A 

10 ' 

X 

10 ' 

06:41:06.051 

+10:18:59.80 

2010Dec 

370 

L142.01+1.77 

A 

10 ' 

X 

10 ' 

03:27:29.910 

+58:43:55.31 

2011Dec 

305 

L151.61-0.24 

A 

10 ' 

X 

10 ' 

04:11:07.205 

+51:09:21.54 

2011Dec 

266 

L154.37+2.58 

A 

10 ' 

X 

10 ' 

04:36:18.153 

+51:11:02.11 

2011Dec 

266 

L169.18-0.89 

A 

10 ' 

X 

10 ' 

05:13:26.754 

+37:27:38.75 

2011Dec 

232 

L172.88+2.27 

A 

10 ' 

X 

10 ' 

05:36:53.085 

+36:10:29.56 

2011Dec 

315 

L211.53-19.27 

A 

10 ' 

X 

10 ' 

05:39:57.873 

-07:27:48.36 

2011Dec 

254 

Llll.42+0.76 

A 

10 ' 

X 

10 ' 

23:12:52.873 

+61:24:31.27 

2012 Sep 

375 

Llll.88+0.99 

A 

10 ' 

X 

10 ' 

23:15:45.489 

+61:47:37.40 

2012 Sep 

226 

L134.20+0.75 

A 

10 ' 

X 

10 ' 

02:28:05.671 

+61:29:25.54 

2012 Sep 

525 

L189.85+0.50 

A 

10 ' 

X 

10 ' 

06:09:19.407 

+20:39:39.18 

2012 Sep 

506 

L189.94+0.34 

A 

10 ' 

X 

10 ' 

06:08:53.761 

+20:30:07.72 

2012 Sep 

590 

L206.60-16.37 

A 

10 ' 

X 

10 ' 

05:41:47.488 

-01:57:59.97 

2012 Sep 

563 


a A: BOX SCAN 571.429 600 40.0 45. B: BOX SCAN 345.600 360 30.0 45. C: BOX SCAN 228.571 240 20.0 45. D: BOX SCAN 606.1 636.4 60.0 
45. E: BOX SCAN 169.7 167.5 60.0 45. F: BOX SCAN 144.6 148.0 60.0 45. G: BOX SCAN 144.6 148.0 60.0 45. H: BOX SCAN 707.1 669.9 60.0 
45. I: BOX SCAN 1060.7 1157.1 60.0 45. The configuration BOX SCAN X Y R A indicates that the map has a size of X(arcsec) X Y(arcsec), a 
scan rate R (arcsec/sec), and a scanning angle A (deg). 

k Representative noise is estimated on the central area of each map, avoiding edges (see Section 4.1). 

C Additional SHARC-II 350 /im continuum data on the Galactic center were presented by Bally et al. (2010). 



Table 3 

Properties of sources recovered in the 350 /i m maps 
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Table 4 

Sources from BGPS catalog with large number of high-resolution 
associated substructures. 


BGPS V2.1 name 

Peak 

emission 

of 350 /im 

substructures 

Total 41 

Qo’r'ms lOcr^TT^^ 2i)(T'r'ms 

BGPSv2_G213.705-12.603 

34 

24 

15 

12 

BGPSv2_G034.256+00.154 

27 

18 

15 

8 

BGPSv2_G133.716+01.220 

25 

15 

9 

4 

BGPSv2_G000.014-00.017 

23 

9 

5 

1 

BGPSv2_G029.916-00.045 

22 

18 

8 

3 

BGPSv2_G029.958-00.017 

19 

13 

6 

4 

BGPSv2_G359.867-00.083 

16 

12 

8 

4 

BGPSv2_G081.477+00.020 

16 

11 

6 

0 

BGPSv2_G359.982-00.069 

14 

12 

7 

1 

BGPSv2_G359.946-00.045 

14 

5 

1 

0 

BGPSv2_G030.751-00.051 

14 

8 

7 

5 

BGPSv2_G081.753+00.593 

13 

9 

7 

5 

BGPSv2_G076.359-00.600 

13 

11 

7 

1 

BGPSv2_G030.786-00.025 

13 

12 

11 

3 

BGPSv2_G203.223+02.076 

12 

8 

1 

0 

BGPSv2_G192.598-00.049 

12 

7 

5 

1 

BGPSv2_G203.320+02.058 

11 

7 

6 

4 

BGPSv2_G023.437-00.183 

11 

5 

2 

1 

BGPSv2_G023.273-00.211 

11 

7 

1 

0 

BGPSv2_G081.721+00.573 

10 

9 

8 

7 

BGPSv2_G030.690-00.043 

10 

4 

3 

0 

BGPSv2_G024.493-00.039 

10 

7 

6 

3 


a Total number of associated substructures (compact and faint) found on 
SHARC-II maps per each BGPS parental clump. See Section 4.5 for details. 
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Table 5 

BGPS V2.1 sources fluxes and fitted dust temperature. 


BGPS source 

1.1 mm 

Convolved 350 /im a 


Temperature 




Flux 40" 

Flux 40" 

/3 = 

1.0 

/3 = 

1.7 

/? = 

2.0 


(Jy) 

(Jy) 

(K) 

(K) 

(K) 

BGPSv2_G024.743+00.179 

0.51T0.11 

17.82i0.64 

>1000.0 

>1000.0 

103.1 

23.5 

31.1 

18.3 

16.7 

19.9 

14.1 

BGPSv2_G024.745+00.161 

0.60T0.12 

18.42i0.66 

351.7 

>1000.0 

58.2 

20.2 

24.8 

16.6 

15.1 

17.3 

13.1 

BGPSv2_G024.757+00.091 

1.36=b0.12 

61.75T2.19 

>1000.0 

>1000.0 

>1000.0 

34.9 

42.9 

29.2 

21.2 

23.6 

19.1 

BGPSv2_G024.759+00.065 

0.43T0.09 

8.47T0.38 

33.3 

52.7 

23.4 

14.2 

16.2 

12.3 

11.6 

12.9 

10.4 

BGPSv2_G024.760+00.163 

0.34T0.13 

4.80T0.18 

20.8 

31.9 

14.2 

11.7 

14.0 

9.5 

10.0 

11.5 

8.4 

BGPSv2_G024.773+00.125 

0.22i0.09 

1.67T0.07 

12.9 

16.2 

9.9 

8.9 

10.3 

7.5 

8.0 

8.9 

6.8 

BGPSv2_G024.791+00.083 

10.07+0.64 

556.65+19.73 

>1000.0 

>1000.0 

>1000.0 

58.6 

78.9 

46.3 

27.1 

30.2 

24.4 

BGPSv2_G024.795+00.101 

3.14T0.22 

142.57T4.87 

>1000.0 

>1000.0 

>1000.0 

34.8 

41.0 

30.1 

21.2 

23.0 

19.5 

BGPSv2_G024.807+00.039 

0.32T0.09 

5.84+0.23 

29.1 

49.6 

19.4 

13.5 

16.0 

11.3 

11.2 

12.7 

9.7 

BGPSv2_G024.815+00.189 

0.10i0.09 

-0.75+0.05 


— 


— 


— 

BGPSv2_G024.817+00.129 

0.48T0.09 

6.69T0.27 

20.7 

25.5 

17.0 

11.6 

12.8 

10.5 

9.9 

10.7 

9.1 

BGPSv2_G024.824+00.181 

0.17A0.09 

1.89T0.07 

16.5 

25.0 

10.7 

10.3 

12.7 

7.9 

9.0 

10.7 

7.2 

BGPSv2_G024.850+00.085 

0.90i0.11 

36.60il.32 

>1000.0 

>1000.0 

>1000.0 

28.9 

35.9 

23.7 

19.0 

21.5 

16.8 

BGPSv2_G024.863+00.145 

0.16T0.09 

1.08zb0.04 

12.2 

16.4 

8.4 

8.6 

10.3 

6.7 

7.7 

9.0 

6.1 

BGPSv2_G025.329-00.196 

0.56T0.12 

3.95T0.16 

12.4 

13.9 

10.9 

8.7 

9.4 

8.0 

7.8 

8.3 

7.2 

BGPSv2_G025.339-00.170 

0.19T0.12 

0.76i0.04 

9.5 

11.7 

7.0 

7.3 

8.4 

5.8 

6.6 

7.5 

5.4 

BGPSv2_G025.355-00.156 

0.19T0.13 

1.49zb0.08 

13.2 

19.5 

8.0 

9.1 

11.3 

6.4 

8.1 

9.7 

5.9 

BGPSv2_G025.355-00.190 

1.46=b0.15 

40.94±1.53 

116.1 

738.0 

61.5 

18.6 

20.6 

16.8 

14.2 

15.3 

13.2 

BGPSv2_G025.382-00.182 

3.23T0.24 

115.44i4.01 

>1000.0 

>1000.0 

>1000.0 

24.0 

26.7 

21.8 

17.0 

18.1 

15.9 

BGPSv2_G025.399-00.140 

4.56A0.31 

152.14i5.34 

>1000.0 

>1000.0 

405.7 

22.1 

24.1 

20.3 

16.1 

17.0 

15.1 

BGPSv2_G025.405-00.256 

0.31T0.15 

3.14i0.11 

15.8 

23.1 

10.6 

10.1 

12.3 

7.9 

8.8 

10.4 

7.1 

BGPSv2_G025.413-00.176 

0.65T0.12 

17.80i0.64 

98.2 

>1000.0 

43.3 

18.2 

21.6 

15.4 

14.0 

15.8 

12.4 

BGPSv2_G025.455-00.210 

1.64=b0.17 

75.41i2.61 

>1000.0 

>1000.0 

>1000.0 

35.9 

46.0 

29.1 

21.5 

24.4 

19.1 

BGPSv2_G025.467-00.126 

0 .12T0.12 

0.33i0.06 

8.4 

11.2 

6.7 

8.2 

6.1 

7.4 

BGPSv2_G025.477-00.136 

0 .10i0.12 

1.73i0.06 

27.0 

>1000.0 

13.1 

25.9 

10.9 

17.8 


Note. — (This table is available in its entirety in a machine-readable form in the online journal. A 
portion is shown here for guidance regarding its form and content.) 

a Fluxes were obtained from SHARC-II maps convolved to match a 33" beam size (see Section 4.6). 
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Flux Density in 40" apertures (Jy) 

Figure 1 . Distribution of 40" aperture flux density of the complete sample of 8594 sources from the BGPS V2.1 catalog (blue line), 
along with those 619 sources contained in the SHARC-II maps (shaded black line). The red line shows the distribution of the 104 sources 
considered as representative targets (see Section 2.1). 
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Figure 2. Error map obtained from CRUSH reduction for L133.71+1.21. The image unit is Jy beam -1 (beam size of 8.5"). The measured 
noise for the L133.71+1.21 map is 0.338 Jy beam -1 . The average representative noise for the complete sample of maps is 478 mjy beam -1 
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Figure 3. Map at 350 obtained with SHARC-II toward L133.71+1.21, corresponding to the W3 Main star forming region. Beam size 
of image is shown in the bottom left corner. Contour levels represent 3cr, 6cr, 10cr ,15cr, 30cr, 50cr and lOOcr, with a rms noise cr = 338 mjy 
beam -1 . 
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Figure 4. Example of the extraction of sources with Bolocat. Left: Map of L030.70-0.07 with SHARC-II at 350 fim (FWHM beam 
of S.S"). Right: The same region mapped with Bolocam at 1.1 mm (FWHM beam of 33"). Some structures of Bolocam have several 
substructures mapped at the better resolution of SHARC-II. Also, the figure shows BGPS sources with no counterpart in the 350 /im map 
(see Section 4.2 for a description of Bolocat parameters used on the extraction of sources). 
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Amplitude [Jy/Beam] 


Figure 5. Completeness fraction derived from fake source tests. The tests were performed in nine SHARC-II maps with different noise 
levels. The plot shows the results of those detections scaled to the mean noise of the catalog maps (0.48 Jy beam -1 ). The fraction of 
sources recovered is plotted as a function of input source intensity. The vertical dotted lines indicate {1, 2, 3, 4, 5, 6}cr. The dashed 
line with diamonds corresponds to fake sources across the whole map, while the continuous line with squares corresponds to sources at 
a distance <3.5' from the center of the map. The results indicate that the catalog is complete at the >99% limit for sources with flux 
densities > 6cr. 



Figure 6. Recovery of source flux density for simulated resolved sources (FHWM = 2.7#6 eam = 23 // ) in SHARC-II maps. Results of the 
calibrated flux recovery are shown in units of Jy. Small differences are found between the 40" aperture and the integrated flux, with an 
under estimation of the integrated flux of ~ 10%. 
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Figure 7. Recovery of source size for simulated observations in SHARC-II maps. The vertical lines at 8.5" and 48" indicate, respectively, 
the beam size of the images, and the largest size recovered by Bolocat for identified objects in our sample of maps. The largest major-axis 
recovered for identified sources is 65". 



Figure 8. Fraction of blended fake objects as a function of source pair separation (solid line). The vertical red line shows the FWHM beam 
size of the 350 p,m maps. Smaller input sources (15"—24" in size) are represented by the dotted line, and large input sources (24"—34" in 
size) are represented by segmented line. 
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Figure 9. Flux distribution for objects identified in the 107 SHARC-II maps. Three methods of flux recovery are represented in the 
image: aperture photometry at 20" (red long dash) and 40" (blue short dash) in diameter toward the peak, and integrated emission of the 
source (solid black line). 



0.35 mm sources radius (") 0.35 mm sources aspect ratio 


Figure 10. Distribution of the deconvolved radii and aspect ratio of sources in the SHARC-II maps. The average and median values of 
the radius distribution are 15" and 14", respectively. The average and median values of the aspect ratios are respectively 1.53 and 1.45. 
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Figure 11. Upper left: Results of Bolocat on the 350 fi m map of L133.71+1.21, corresponding to the W3 Main region. The different colors 
indicate the peak signal-to-noise of each source (above 10cr in red, below that limit in blue). Sources in red are associated with compact 
and dense core-type structures, the blue sources are related to low emission, field material or possible filamentary structure. Upper right: 
Grey color scale of the 1.1 mm toward the same region. Black regions show the result of Bolocat on the 1.1 mm map, and blue and red 
regions the results of Bolocat on the 350 [im map. The central V2.1 source mask shows large number of substructures (both strong and 
faint). Lower panels are similar to the upper panels for L081.48 + 0.00, in the Cygnus X region. The comparison of source masks reveal a 
large population of BGPS V2.1 sources located inside the 350 /im maps, but not associated with any source at high-resolution. Details on 
the correlation between 350 /am sources and their parental BGPS clump are described in Section 4.5. 
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Figure 12. Distribution of 350 jtmi sources per 1.1 mm source from BGPS. The top panel shows the cumulative number. Half of the 1.1 
mm parental clumps are associated with one-or-two substructures identified at high-resolution. 
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Figure 13. Upper: Comparison between integrated flux in 40" apertures obtained at 350 /im on SHARC-II and Herschel images toward 
l = 30°. The points represent the 213 high-resolution sources recovered by Bolocat in SHARC-II maps, and the colors show different limits 
of peak signal-to-noise on those sources. Bottom: Comparison between integrated flux in 40" apertures for 350 /im images from SHARC-II 
and Herschel image on a sample of 102 BGPS sources located toward l = 30°. Both SHARC-II and Herschel were convolved to match the 
resolution of the Bolocam maps (33"). Unfilled circles represent BGPS sources with no associated source detected at 350 /im, and filled 
circles show the 1.1 mm parental clumps with only faint objects (black) and with compact substructures (red). The error bars in the flux 
ratio considered 20% uncertainties in F^erscheU and 30% in Fsharc • 
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Right Ascension (J2000) 



Figure 14. Example of the 350 fi m continuum emission from SHARC-II (upper left) and Herschel /SPIRE (upper right) toward the region 
l = 28.95°, b = —0.05°. Beam sizes are shown in the bottom left corner of each map. The contours of the SHARC-II image correspond 
to 3cr, 6cr, 10(7, 15cr, 30cr, 50cr and 100cr, with an rms noise a — 514 mJy beam -1 . Contours of the Herschel image start from 10cr, with 
increasing steps of 5cr (cr = 164 MJy sr -1 ). The smaller beam size of SHARC-II images allow us to recover structures not identified by 
Herschel with lower resolution. Bottom: SHARC-II image convolved to match the 24.9 /x FWHM beam size of Herschel image. Contours 
represent 1, 3, 6, 10, 15, 30 and 50 times the rms noise of the image. The lcr contour level resembles the emission observed at the 15(7 
white contour level of the Herschel image. The 15cr level in the Herschel map sets then a limit of the diffuse emission filtered out in the 
SHARC-II maps (see Section 4.7). 
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Figure 15. Distribution of the flux ratio between the emission at 0.35 mm and 1.1 mm for our sample of 574 BGPS V2.1 sources with 
reliable flux values. Nearly 94% of the sources have a flux ratio lower than 70.2 
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Figure 16. Models of dust temperature as a function of the flux ratio between continuum emission at 350 [im and 1.1 mm, in the optically 
thin limit. The three models correspond to different opacity laws, (3 = 1.0, (3 = 1.7 and (3 = 2.0. The error bars for the (3 = 1.7 model show 
the expected uncertainty in temperature determination when the flux ratio has an error of 10%. 
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Figure 17. Distribution of temperatures of the BGPS sources with good fit in the determination of temperature, using the three different 
opacity models: spectral index (3=1.0 (blue line), (3=2.0 (green dashed line), and (3=1.7 (shaded black line). The vertical red dashed line 
at T=13.3 K, represents the median value for the distribution with (3=1.7. 
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Figure 18. Temperature determined for the BGPS sample of sources as a function of their 40" aperture flux, for a spectral index (3 = 1.7. 
Only sources with good fits in the determination of temperature (T< 50 K) are shown in this figure. There are few BGPS sources with 
high flux (Fi.i mm > 4 Jy) and low temperature (T< 20 K); a description of them is presented in Section 5.1. 
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Figure 19. Upper: Comparison between estimated temperatures with 350 /im and 1.1 mm continuum maps, TsharCi and temperatures 
obtained from SED fitting by Battersby et al. (2011), T^jerschel: as a function of 40" aperture flux density at 1.1 mm for a sample of 91 BGPS 
sources toward l = 30°. The colors shown in the legend represent different ranges for Tsharc • The mean value of the ratio, weighted by 
errors, is 0.50T0.01 for sources with Tsharc < 50 K. Bottom: Simulated recovery of sources immersed in a uniform background emitting 
as a modified blackbody with temperature Tj^=25 K, and spectral index = 1.7. The input sources have temperatures and spectral 
index in the ranges 5 < T source < 50 K, and 1.5 < /3 source < 2.5, which are characteristic of clump-like features. The temperature ratio 
Tsource/T source + b g of simulated sources immersed in background reproduce qualitatively the observed ratio TsHARc/Tnerschel of the 
upper panel. 



Ratio Tgas / Tdust 

Figure 20. Distribution of the ratio of gas temperature determined from NH 3 observations Tgas: and fitted dust temperatures for 

the sample of 49 BGPS sources with Tgas < 30 K and Td ust < 50 K, ranges in which both values are relatively well determined. The 
average and median values of this distribution are 0.88 and 0.76, respectively. 
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Temperature (K) 



Figure 21. Peak main-beam temperature, as a function of dust temperature, for the sample of BGPS sources from the spectroscopy 
catalog by Shirley et al. (2013) that are contained in the SHARC-II maps. The number of sources with HCO + emission is 250, and the 
number of sources with N 2 H+ emission is 199. 



Figure 22. Ratio between the mass estimated using color temperatures determined from 350 /im and 1.1 mm images (T c i urnp ), and the 
mass estimated using a single fiducial temperature ( Tf ix ). For the sample of 514 V2.1 BGPS sources with good fits in the determination 
of temperature, using Tf^ x = 20 K, the average value of the source masses will be underestimated by a factor of ~2.4. 
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Figure 23. Comparison between the total estimated mass from 350 jttm substructures in each clump, and the mass estimated at 1.1 mm 
for that parental BGPS clump, as a function of dust temperature (left) and the integrated flux of the BGPS sources (right). Black points 
consider for M su b structures the flux contribution from all high-resolution sources in the parental clumps, and have a weighted averaged 
mass fraction ^(^-substructures)/Mciump of 0.22=1=0.01. Red points only consider the flux contribution from strong, compact substructures, 
and the averaged mass fraction, weighted by errors, of this distribution is 0.19=b0.01. 



0.001 0.010 0.100 1.000 10.000 
Surface density parental clump (g cm" 2 ) 


Figure 24. Estimation of surface density between substructures and parental clumps. Solid line represent the locus of equality between the 
structures. Black crosses show the bulk of the sources, with temperatures between 10—50 K. Blue marks show sources with low temperature 
(T< 10 K), likely to overestimating the surface density. Grey marks show sources with high temperature (T> 50 K). Red points show 
only those compact, core-type substructures found in the parental clumps, with detections above 10cr. Most of these compact sources 
have higher surface 2-3 times larger than their parental clump. Several of these compact sources have surface densities above theoretical 
thresholds of 0.2 g cm -2 (Butler & Tan, 2012) or 1.0 g cm -2 (Krumholz &; McKee, 2008) for massive star formation (see Section 6.1). 
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APPENDIX 

CALIBRATORS 

Table 6 presents in detail the flux calibration and the observed values of r at 225 GHz for each observation run. Flux 
conversion factors for one beam ( C^eam ), 20" and 40" aperture conversion factors on SHARC-II maps (C 20 and C 40 ), 
and 40" aperture conversion factors for SHARC-II maps convolved to a resolution of 33" are shown in the table. 


Table 6 

Calibrators 


Date 

Calibrator 

T225 GHz 

Ffoeam 
(Jy beam -1 

c 20 

(Jy vV- 1 ) 

c 40 

(Jy MV 1 ) 

P 33 " 

U 40 1 

(Jy mV- 1 ) 

2006 Jun 11 

Mars 

0.071 

10.86 

0.39 

0.31 

0.46 


Mars 

0.071 

10.15 

0.42 

0.30 

0.46 


Mars 

0.071 

11.63 

0.40 

0.33 

0.49 

2006 Jun 23 

Uranus 

0.061 

8.55 

0.28 

0.24 

0.36 


Uranus 

0.062 

8.02 

0.28 

0.23 

0.34 

2007 Jul 9 

Mars 

0.068 

7.56 

0.27 

0.22 

0.33 


Uranus 

0.058 

6.55 

0.22 

0.18 

0.27 


Uranus 

0.068 

7.57 

0.25 

0.21 

0.31 


Neptune 

0.058 

6.76 

0.22 

0.19 

0.28 


Neptune 

0.064 

6.42 

0.21 

0.18 

0.27 


Neptune 

0.060 

7.18 

0.23 

0.20 

0.30 

2007 Oct 24 

Uranus 

0.053 

7.80 

0.26 

0.22 

0.32 


Uranus 

0.053 

8.02 

0.27 

0.23 

0.34 


Neptune 

0.055 

9.29 

0.31 

0.26 

0.39 


Neptune 

0.055 

9.20 

0.30 

0.26 

0.38 

2007 Oct 27 

Uranus 

0.052 

7.18 

0.23 

0.20 

0.30 


Uranus 

0.047 

7.80 

0.25 

0.22 

0.33 


Uranus 

0.047 

7.29 

0.24 

0.20 

0.30 

2009 Sep 10 

Uranus 

0.055 

5.80 

0.19 

0.16 

0.24 


Mars 

0.062 

6.31 

0.23 

0.18 

0.27 

2009 Dec 13 

Uranus 

0.050 

6.34 

0.22 

0.18 

0.26 


Uranus 

0.050 

6.13 

0.23 

0.18 

0.26 


Mars 

0.065 

5.73 

0.20 

0.17 

0.24 


Mars 

0.065 

5.67 

0.20 

0.16 

0.24 

2009 Dec 14 

Uranus 

0.068 

6.01 

0.22 

0.17 

0.25 


Mars 

0.058 

5.28 

0.22 

0.16 

0.24 

2009 Dec 29 

Uranus 

0.060 

5.84 

0.20 

0.16 

0.24 


Uranus 

0.052 

5.85 

0.19 

0.16 

0.24 


Mars 

0.062 

5.72 

0.22 

0.17 

0.26 


Mars 

0.062 

5.78 

0.22 

0.17 

0.26 


Mars 

0.056 

5.80 

0.22 

0.17 

0.26 

2009 Dec 30 

Uranus 

0.050 

7.06 

0.24 

0.19 

0.28 


Uranus 

0.051 

7.08 

0.24 

0.19 

0.28 


Mars 

0.047 

5.70 

0.22 

0.17 

0.25 


Mars 

0.052 

5.68 

0.22 

0.17 

0.25 


Mars 

0.040 

5.66 

0.22 

0.17 

0.25 


Mars 

0.040 

5.81 

0.24 

0.18 

0.26 

2009 Dec 31 

Uranus 

0.051 

6.43 

0.22 

0.18 

0.27 


Uranus 

0.049 

6.32 

0.22 

0.18 

0.26 


Uranus 

0.050 

6.26 

0.22 

0.18 

0.27 


Mars 

0.048 

4.87 

0.19 

0.15 

0.22 


Mars 

0.048 

4.79 

0.20 

0.15 

0.22 

2010 Jan 01 

Mars 

0.045 

5.69 

0.22 

0.17 

0.25 

2010 Jul 24 

Neptune 

0.063 

19.30 

0.61 

0.54 

0.80 


Uranus 

0.070 

26.18 

0.87 

0.72 

1.07 


Uranus 

0.075 

18.00 

0.76 

0.52 

0.81 


Uranus 

0.056 

11.68 

0.58 

0.36 

0.55 


Uranus 

0.065 

13.98 

0.55 

0.39 

0.61 

2010 Jul 25 

Neptune 

0.077 

9.63 

0.32 

0.27 

0.41 


Uranus 

0.042 

6.09 

0.20 

0.17 

0.25 


Uranus 

0.039 

5.79 

0.19 

0.16 

0.24 


Uranus 

0.046 

7.12 

0.25 

0.20 

0.30 

2010 Jul 28 

Neptune 

0.056 

6.76 

0.21 

0.19 

0.28 


Uranus 

0.056 

6.57 

0.25 

0.18 

0.27 


Uranus 

0.046 

6.37 

0.20 

0.18 

0.26 

2010 Jul 31 

Neptune 

0.050 

8.58 

0.28 

0.23 

0.35 


Neptune 

0.050 

7.71 

0.26 

0.21 

0.32 


Neptune 

0.050 

7.82 

0.27 

0.22 

0.33 


Neptune 

0.050 

8.22 

0.27 

0.23 

0.34 


Uranus 

0.058 

6.67 

0.23 

0.19 

0.28 

2010 Aug 01 

Neptune 

0.036 

7.42 

0.24 

0.20 

0.30 


Neptune 

0.060 

7.67 

0.25 

0.21 

0.31 


Neptune 

0.050 

7.46 

0.25 

0.20 

0.30 
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Table 6 — Continued 


Date 

Calibrator 

r 225 GHz 

C beam 

(Jy beam -1 iV -1 ) 

C20 

(Jy uV- 1 ) 

c 40 

(Jy uV - 1 ) 

P 33" 

U 40 1 

(Jy mV- 1 ) 


Uranus 

0.053 

7.06 

0.24 

0.19 

0.29 

2010 Dec 05 

Uranus 

0.060 

5.75 

0.20 

0.17 

0.24 


Uranus 

0.030 

5.96 

0.21 

0.17 

0.25 


Uranus 

0.030 

6.66 

0.23 

0.19 

0.27 


Uranus 

0.040 

6.38 

0.21 

0.18 

0.26 

2010 Dec 06 

Uranus 

0.030 

7.07 

0.25 

0.19 

0.29 


Uranus 

0.050 

7.07 

0.26 

0.20 

0.30 


Uranus 

0.038 

7.26 

0.24 

0.19 

0.30 


Uranus 

0.038 

7.04 

0.25 

0.19 

0.29 


Uranus 

0.068 

6.78 

0.26 

0.19 

0.29 


Uranus 

0.049 

6.93 

0.23 

0.19 

0.28 

2011 Dec 23 

Uranus 

0.036 

11.35 

0.38 

0.32 

0.49 


Uranus 

0.071 

10.06 

0.37 

0.29 

0.43 


Uranus 

0.070 

11.06 

0.38 

0.31 

0.45 

2012 Sep 14 

Uranus 

0.050 

12.73 

0.42 

0.35 

0.52 


Uranus 

0.057 

12.88 

0.44 

0.36 

0.53 


Uranus 

0.057 

12.60 

0.42 

0.35 

0.52 


Uranus 

0.062 

17.32 

0.57 

0.47 

0.70 


Uranus 

0.065 

17.37 

0.59 

0.48 

0.72 

2012 Sep 21 

Uranus 

0.070 

17.47 

0.56 

0.49 

0.73 


Uranus 

0.094 

17.79 
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SHARC-II MAPS 

We present in Figure 25 our complete sample of 107 regions mapped with SHARC-II. For each of them, the image 
shows grey scale maps with emission contours at 3a, 6a, lOcr, 15a, 30a, 50a, and 100a, with a the rms noise. The value 
of a in units of mJy beam -1 is stated below each map. The image also shows the high-resolution sources extracted 
with Bolocat. Red contours represent sources with signal-to-noise values above 10 (compact sources), and blue contours 
sources below this limit (faint sources). 
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(25.1) L359.85+0.00 map, a r ms = 994 mjy beam 1 . 
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(25.2) L000.00+0.00 map, a r ms = 1453 mJy beam 
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(25.3) L000.15+0.00 map, a r ms = 1438 mjy beam 1 . 


Figure 25. Continuum maps at 350 fi m of the 107 regions included in our catalog. The beam size of the maps is . Left: Contours 
are drawn at 3a, 6a, 10 a, 15 a, 30 a, 50 a, and 100cr. Noise level a r ms is indicated in the figure for each map. Right: Recovery of substruc¬ 
tures with Bolocat algorithm. Red and blue regions represent the position of dense, compact cores-type sources and faint substructures, 
respectively. 
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(25.4) L000.30+0.00 map, cr rms = 1086 mjy beam 1 . 
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(25.5) L029.95-0.05 map, cr rrns = 514 mjy beam 1 . 
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(25.6) L030.00+0.00 map, <J r ms = 462 mjy beam 1 . 


39 


Figure 25. Continuation 
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(25.7) L030.15+0.00 map, cr rms = 426 mJy beam 1 . 
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(25.8) L030.30+0.00 map, cr rms = 496 mJy beam -1 . 
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(25.9) L030.45+0.00 map, cr rms = 603 mJy beam 


Figure 25. Continuation 
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(25.11) L030.70-0.07 map, cr rms = 510 mjy beam 1 . 
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(25.12) L030.80-0.05 map, cr rms = 529 mjy beam 1 . 


Figure 25. Continuation 
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(25.13) L030.88+0.13 map, cr r ms = 367 mJy beam 1 . 
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(25.14) L034.26+0.15 map, cr rms = 395 mJy beam 1 . Additional contour is drawn at 200cr and 
400cr. 
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(25.15) L076.16+0.09 map, cr r ms = 249 mJy beam 1 . 


Figure 25. Continuation 
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(25.16) L077.93+0.02 map, cr r ms = 389 mJy beam 1 . 
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(25.17) L078.14-0.32 map, cr rms = 408 mJy beam 1 . 
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Figure 25. Continuation 
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(25.19) L079.28+0.30 map, cr rms = 558 mJy beam 1 . 
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(25.20) L080.92-0.ll map, a r ms = 724 mJy beam 1 . 
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(25.21) L081.45+0.04 map, a r ms = 274 mJy beam 1 . 


Figure 25. Continuation 
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(25.22) L081.68+0.54 map, cr r ms = 802 mJy beam 
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(25.23) L082.55+0.14 map, cr r ms = 255 mJy beam 
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Figure 25. Continuation 
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(25.25) L076.12-0.24 map, 



(7rms = 388 mjy beam 1 . 
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Figure 25. Continuation 








The Bolocam Galactic Plane Survey. XI 


47 


L078.17-0.31 



Right Ascension (J2000) 


(25.28) L078.17-0.31 map, 


L078.17-0.31 



o'rms = 800 mjy beam 1 . 
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(25.29) L079.62+0.49 map, cr rms = 1489 mjy beam 1 . 
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(25.30) L079.11-0.35 map, a rrns = 1226 mjy beam - 1 . 


Figure 25. Continuation 
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(25.31) L080.86+0.38 map, cr r ms = 1874 mjy beam 1 . 




(25.32) L110.11+0.05 map, a rms = 288 mjy beam 1 . 
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(25.33) Llll.62+0.38 map, cr rms = 342 mjy beam - 1 . 


Figure 25. Continuation 
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(25.34) L134.28+0.86 map, cr r ms = 275 mJy beam 1 . 
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(25.36) L136.38+2.27 map, cr rms = 257 mJy beam 1 . 


Figure 25. Continuation 
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(25.37) L136.83+1.07 map, cr rms = 250 mjy beam -1 . 
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(25.38) Llll.28-0.66 map, cr rms = 286 mjy beam 1 . 
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(25.39) Llll.54+0.78 map, cr r ms = 496 mjy beam 1 . 


Figure 25. Continuation 
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(25.40) L133.71+1.21 map, cr rms = 338 mJy beam 
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(25.41) L133.95+1.06 map, a rrns = 338 mJy beam -1 . Additional contour is drawn at 200cr. 


Llll.26-0.77 


Llll.26-0.77 




(25.42) Llll.26-0.77 map, 


= 538 mJy beam 
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(25.45) Llll.88+0.82 map, cr r ms = 473 mjy beam 1 . 


Figure 25. Continuation 
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(25.49) L137.69+1.46 map, cr r ms = 667 mjy beam -1 . 
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(25.50) L138.30+1.56 map, cr r ms = 205 mjy beam 1 . 
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(25.51) L138.48+1.63 map, cr r ms = 266 mjy beam 


Figure 25. Continuation 
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(25.53) L173.17+2.35 map, cr r ms = 320 mJy beam 1 . 
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(25.54) L173.47+2.43 map, a r ms = 487 mJy beam 
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(25.55) L173.57+2.44 map, cr rms = 421 mJy beam 1 . 
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(25.56) L173.62+2.81 map, cr rms = 359 mJy beam 1 . 
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(25.58) L173.76+2.67 map, cr r ms = 301 mJy beam 1 . 
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(25.59) L188.79+1.03 map, a rms = 259 mJy beam 1 . 
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(25.60) L188.95+0.88 map, a rms = 268 mJy beam 1 . 
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(25.64) L189.85+0.39 map, a rms = 505 mjy beam 1 . 
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(25.67) L192.60-0.05 map, cr r ms = 271 mjy beam 1 . 
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(25.71) L196.42-l.66 map, cr rms = 470 mJy beam 1 . 
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(25.72) L203.23+2.06 map, cr r ms = 269 mJy beam 1 . 
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(25.73) L203.35+2.03 map, cr r ms = 314 mJy beam 
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(25.74) L213.71-12.62 map, cr rms = 421 mJy beam 
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(25.75) L217.37-0.07 map, cr rms — 704 mJy beam 1 . 
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Figure 25. 
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(25.82) L024.78+0.12 map, cr r ms = 363 mJy beam 




(25.83) L025.40-0.18 map, cr rms = 482 mJy beam 1 . 
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(25.84) L030.61+0.16 map, cr r ms = 349 mJy beam - 1 . 
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(25.86) L081.ll-0.16 map, a r ms = 341 mJy beam -1 . 
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(25.87) L081.28+1.01 map, a r ms = 250 mJy beam 
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(25.88) L081.39+0.73 map, cr r ms = 341 mjy beam 1 . 
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(25.89) L081.48+0.00 map, cr rms = 353 mJy beam 1 . 
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(25.90) L081.76+0.60 map, cr rms = 329 mjy beam 1 . Additional contour is drawn at 200cr. 
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(25.91) L031.41+0.31 map, a r ms = 791 mJy beam 1 . 
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(25.92) L081.88+0.77 map, a rms = 717 mJy beam -1 . 
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(25.93) L183.40-0.58 map, a r ms = 255 mJy beam 
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(25.94) L189.78+0.33 map, cr r ms = 294 mJy beam 
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(25.95) L202.58+2.42 map, a r ms = 370 mJy beam 
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(25.96) L142.01+1.77 map, a r ms = 305 mJy beam 


69 
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(25.97) L151.61-0.24 map, cr r ms = 266 mjy beam 
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(25.98) L154.37+2.58 map, a r ms = 266 mJy beam 
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(25.99) L169.18-0.89 map, cr rms = 232 mjy beam 1 . 
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(25.100) L172.88+2.27 map, a r ms = 315 mJy beam -1 . 
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(25.101) L211.53-19.27 map, a rrns = 254 mJy beam 1 . 
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(25.102) Llll.42+0.76 map, a rrns = 375 mJy beam 1 . 
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(25.103) Llll.88+0.99 map, arms = 226 mjy beam 1 . 
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(25.104) L134.20+0.75 map, a r ms = 525 mjy beam 1 . 
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(25.105) L189.85+0.50 map, cr rms = 506 mjy beam 1 . 
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(25.106) L189.94+0.34 map, cr rms = 590 mjy beam 
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(25.107) L206.60-16.37 map, cr rms = 563 mjy beam 


73 


Figure 25. Continuation 






74 


Merello et al. 


HI-GAL SPIRE IMAGES AT 350 /im FOR l = 30° 

In this section we indicate the locations of the high-resolution sources identified in the 350 fi m SHARC-II maps 
(8.5"beam size) toward l = 30°, and their comparison with structures observed at lower resolution (24.9") in the 
corresponding Herschel/SPIKE maps at the same wavelength and BGPS maps at 1.1 mm at a resolution of 33". 
Figure 26 shows 1.1 mm continuum gray scale maps with green contours representing sources from the BGPS catalog, 
and Herschel grey scale maps with contours representing the emission at lOcr, and increasing steps of 5a, with a=164 
MJy sr -1 . The thick black line on each map represent the area covered by the SHARC-II maps, with the name of those 
maps indicated for each region. Peak positions of detected high-resolution sources from SHARC-II maps are overlaid 
in red and blue crosses, representing sources with signal-to-noise above and below a value of 10, respectively. Most 
sources with high signal-to-noise are well associated with strong emission in the Herschel maps, but several of the low 
signal-to-noise sources are not individually identified and they are part of extended background emission detected in 
Herschel maps but not recovered in SHARC-II maps. 
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(26.1) Thick black contour in Herschel map represents the area covered by maps L029.95-0.05 and L030.00+00. 
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L030.45+0.00. 
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Figure 26. Left: BGPS map at 1.1 mm toward l = 30°. The ellipses show the recovered sources from the BGPS catalog. Right: 
Corresponding 350 p,m maps from Herschel SPIRE. Contours start from 10cr, with increasing steps of 5cr (a =164 MJy sr -1 ) The crosses 
show the peak position of the high-resolution sources recovered in SHARC-II map. Red and blue crosses show substructures with peak 
signal-to-noise above and below 10, respectively. The area covered by SHARC-II maps, indicated under each figure, is represented by the 
thick black contour. 
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(26.3) Thick black contour in Herschel map represents the area covered by maps L030.60+0.00, L030.70-0.07 and L030.80-0.05 
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(26.4) Thick black contour in Herschel map represents the area covered by map L030.88+0.13. 


BGPS 1.1 mm 



Herschel 350 /im 



(26.5) Thick black contour in Herschel map represents the area covered by map L030.61+0.16. 


Figure 26. Continuation. 


















